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Abstract. The Local Cluster Substructure Survey (LoCuSS, Smith et al.) is a systematic multi-wavelength survey of more 
than 100 X-ray luminous galaxy clusters in the redshift range 0.14-0.3 selected from the ROSAT All Sky Survey. We used 
data on 37 LoCuSS clusters from the XMM-Newton archive to investigate the global scaling relations of galaxy clusters. The 
scaling relations based solely on the X-ray data (S-T, S-Y x , P-Yx, M-T, M-Yx, M-M gas , M m -T, L-T, L-Y x , and L-M) obey 
empirical self-similarity and reveal no additional evolution beyond the large-scale structure growth. They also reveal up to 17 
per cent segregation between all 37 clusters and non-cool core clusters. Weak lensing mass measurements are also available in 
the literature for 19 of the clusters with XMM-Newton data. The average of the weak lensing mass to X-ray based mass ratio is 
1.09 + 0.08, setting the limit of the non-thermal pressure support to 9 + 8 per cent. The mean of the weak lensing mass to X-ray 
based mass ratio of these clusters is — 1, indicating good agreement between X-ray and weak lensing masses for most clusters, 
although with 31-51 per cent scatter. The scatter in the mass-observable relations (M-Yx, M-M g2:s , and M-T) is smaller using 
X-ray based masses than using weak lensing masses by a factor of 2. With the scaled radius defined by the K x profile - rj *' x , 
r 50o' Wl ' an d r 5oo i ' we obtain lower scatter in the weak lensing mass based mass-observable relations, which means the origin of 
the scatter is M wl and M x instead of Yx- The normalization of the M-Yx relation using X-ray mass estimates is lower than the 
one from simulations by up to 18-24 per cent at 3cr significance. This agrees with the M-Yx relation based on weak lensing 
masses, the normalization of the latter being ~ 20 per cent lower than the one from simulations at ~ 2cr significance. This 
difference between observations and simulations is also indicated in the M-M gas and M-T relations. Despite the large scatter 
in the comparison of X-ray to lensing, the agreement between these two completely independent observational methods is an 
important step towards controlling astrophysical and measurement systematics in cosmological scaling relations. 
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1. Introduction 

The gravitational growth of fluctuations in the matter density 
distribution can be traced by the evolution of the galaxy clus- 
ter mass function (e.g. Schuecker et al. 2003). The mass func- 
tion of luminous galaxy clusters probes the cosmic evolution 
of large-scale structure (LSS) and is an extremely effective test 
of cosmological models. It is sensitive to the matter density, 
Q m , and the amplitude of the cosmic power spectrum on clus- 
ter scales, erg (e.g. Schuecker et al. 2003). 
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* This work is based on observations made with the XMM-Newton, 
an ESA science mission with instruments and contributions directly 
funded by ESA member states and the USA (NASA). 



The construction of the mass function of galaxy clusters 
for large cosmological cluster samples is based on robust cali- 
bration of cluster mass-observable scaling relations for repre- 
sentative galaxy cluster samples using reliable mass measure- 
ments (e.g. Reiprich & Bohringer 2002; Voit 2005; Stanek et 
al. 2006). 

Deep X-ray observations using XMM-Newton and Chandra 
can recover precisely the intra-cluster medium (ICM) density 
and temperature distributions. XMM-Newton (also Chandra) 
allows us for the first time to measure mass distributions for 
relaxed nearby clusters with statistical uncertainties below 15 
per cent up to rsoo, which yields a tight M500-T relation with 
only 8 per cent scatter (Arnaud et al. 2005). Using Chandra 
data, Vikhlinin et al. (2006a) measured precise mass distribu- 



2 



Zhang et al.: LoCuSS: comparison of observed X-ray and lensing galaxy cluster scaling relations with simulations 



tions up to radii of r25oo, which provides a tight M2500-T re- 
lation for relaxed nearby clusters with its scatter also within 8 
per cent. However, relaxed clusters only represent a small frac- 
tion of the cluster population in the local Universe and an even 
smaller fraction at high redshifts (Jeltema et al. 2005; Smith et 
al. 2005; Vikhlinin et al. 2006b) and the cluster mass estimates 
from X-ray observations are limited by additional physical pro- 
cesses in the ICM. Recent simulations (e.g. Rasia et al. 2004, 
2006; Nagai et al. 2007a) show that the X-ray total mass mea- 
sured via hydrostatic equilibrium can be biased low by up to 20 
per cent for nearby relaxed clusters and therefore the normal- 
ization of the mass-temperature relation can be biased low by 
up to 30 per cent. The bias is more significant (up to ~30 per 
cent) for un-relaxed clusters. 

Gravitational lensing provides a direct probe of the total 
cluster mass without needing to invoke the symmetry and equi- 
librium assumptions required in the X-ray analysis. However 
lensing yields direct constraints on the mass projected along 
the line-of-sight through the cluster and may therefore be prone 
to projection effects. It is therefore very powerful to use both 
X-ray and lensing methods to measure cluster mass and thus 
to investigate systematic uncertainties in mass measurements 
and their impact on cosmological parameter estimation. Early 
X-ray versus strong lensing mass comparisons revealed large 
discrepancies between the methods (Miralda-Escude & Babul 
1995), which were subsequently attributed to the absence of 
a cool core in, and thus presumed more disturbed dynamical 
state of, the discrepant clusters (Allen 1998; Wu 2000). More 
recent analysis of Chandra and Hubble Space Telescope data 
suggest a picture that is broadly consistent with Allen (1998) in 
that disturbed, likely merging, clusters dominate the scatter in 
the mass-temperature relation for cluster cores (R<250 kpc/h) 
when using lensing to infer cluster mass (Smith et al. 2003, 
2005). Wide-field weak lensing studies of cluster scaling re- 
lations have to date generally lacked the statistical precision 
to investigate systematic differences between disturbed and 
undisturbed clusters (Bardeau et al., 2007; Pedersen & Dahle 
2007; Zhang et al., 2007; Hoekstra 2007). However extending 
the mass range of lensing-based scaling relation studies down 
to groups has been shown to help overcome the large scatter 
when trying to measure the slope of, for example, the mass- 
temperature relation (Berge et al. 2007). 

In summary, credible cluster cosmology experiments re- 
quire well calibrated measurements of the shape, scatter and 
evolution of the mass-observable scaling relations based on 
large statistical samples of clusters that are unbiased with re- 
spect to cluster morphology. Elimination of systematic uncer- 
tainties from this calibration demands that the cluster mass 
measurements are cross-checked between independent mass 
measurement techniques, including X-ray and lensing meth- 
ods. The challenges for precision cluster cosmology on the 
per cent level are therefore as follows: (1) construction of 
large (n c ius>100) representative cluster samples, and (2) de- 
tailed high quality X-ray and lensing mass measurements of 
that sample as a basis for calibrating the mass-observable scal- 
ing relations. Despite recent progress on weak lensing searches 
for galaxy clusters (Wittman et al., 2006; Massey et al., 2007; 
Miyazaki et al., 2007), the most efficient method to construct 



a homogeneous sample of massive clusters is to use all sky 
X-ray surveys. The main attractions of this approach are the 
well-defined sensitivity limit of such surveys and the minimal 
biases towards different cluster morphologies. 

The Local Cluster Substructure Survey (LoCuSS, PI: G. P. 
Smith; Smith et al. in prep.) is a systematic multi-wavelength 
survey of galaxy clusters designed to tackle the probems out- 
lined above. LoCuSS provides a morphology-unbiased sample 
of ~100 X-ray luminous galaxy clusters in the redshift range 
0.14<z<0.3 selected from the ROSAT All Sky Survey catalogs 
(Ebeling et al. 1998, 2000; Bohringer et al. 2004). Studies of 
the LoCuSS sample are in progress on (1) gravitational lens 
modeling of the mass distribution in the cluster cores, using 
data from our ongoing Hubble Space Telescope (HST) ob- 
serving program (SNAP:10881, GO:11312; Hamilton-Morris 
et al., in prep.); (2) the scaling relation between cluster mass 
inferred from lensing and the Sunyaev-Zeldovich (SZ) effect Y- 
parameter (Marrone et al., in prep.); (3) the distribution of the 
luminosity gap statistic in 1O 15 M clusters (Khosroshahi et al., 
in prep.); (4) wide-field weak-lensing analysis of clusters ob- 
served with Suprime-CAM on Subaru (Okabe et al., in prep.). 
A wide variety of follow-up observations are also underway 
or imminent, including recent awards of 300 ks on Chandra 
(PID:9800372), 129 ks at 24 fim with Spitzer (PID:40872), 
60 ks on GALEX (Cycle 4, Proposal # 90), and a 0.5 Msec 
Key Program on ESA's Herschel. These space-based obser- 
vations are supported by ongoing ground-based programs at 
Keck, VLT, Subaru, Palomar, NOAO and UKIRT. 

The goal of this paper is to compare X-ray and lensing clus- 
ter mass measurements of LoCuSS clusters for which archival 
X-ray data are available from XMM-Newton, and weak lensing 
masses are available in the literature. As at the end of 2006, 37 
LoCuSS clusters had sufficiently high quality archival XMM- 
Newton data to derive precise temperature distributions for sta- 
tistically reliable X-ray mass estimates. Weak lensing mass 
measurements are also available in the literature for 19 of these 
clusters (including 7 duplicates): 10 clusters from CFH12k data 
(Bardeau et al., 2005; 2007); 15 clusters from observations at 
the Nordic Optical Telescope and UH 88in (Dahle 2006); one 
cluster, Abell 1689, has also been studied using both CFHT 
and Subaru data (Limousin et al., 2007; Halkola et al. 2006; 
Broadhurst et al., 2005). Weak lensing measurements for an en- 
larged sample of LoCuSS clusters are coming, based on high 
quality HST and Subaru data (Okabe et al., in prep.). 

This paper is organised as follows. The XMM-Newton data 
quality and corresponding optimized reduction method are de- 
scribed in § [2j § [3] shows an empirical similarity of the X-ray 
scaling relations for our sample and the difference in the nor- 
malization between mass-observable relations from X-ray ob- 
servations and simulations. In § |4j we used a new method to 
determine the radial scale (jVqq ) to study the scaling relations, 
in which the scatter can be reduced. In § [5] we present X-ray 
versus weak lensing mass calibration using different mass def- 
initions to look for the one which gives least scatter in the 
X-ray versus weak lensing mass comparison. § [6] describes 
the weak lensing mass based mass-observable relations and 
the difference in the normalization between the weak lensing 
mass based mass-observable relations for 19 LoCuSS clus- 
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ters and the relations from simulations. In §|7l we discuss the 
scatter and bias in weak lensing mass estimates, the cluster 
mass comparison between observations and simulations, and 
the variation of the amplitude of cosmic power spectrum with 
the mass calibration. We draw our conclusions in § [8] Unless 
explicitly stated otherwise, we adopt a flat ACDM cosmology 
with the density parameter Q m = 0.3 and the Hubble constant 
Hq = 70 km s Mpc~'. We adopt the solar abundance table 
of Anders & Grevesse (1989). Confidence intervals correspond 
to the 68 per cent confidence level. We apply the Orthogonal 
Distance Regression package (ODRPACK 2.010, e.g. Boggs et 
al. 1987) taking into account measurement errors on both vari- 
ables to determine the parameters and their errors of the fitting 
throughout this paper. We use Monte Carlo simulations for the 
uncertainty propagation on all quantities of interest. 

2. Data description and X-ray mass modeling 

Achieving the goals of this paper (§ [TJ, specifically identify- 
ing systematic differences between mass measurement meth- 
ods, rests on being able to measure statistically reliable cluster 
masses using X-ray data. For this, a precise temperature dis- 
tribution measurement is required. In order to obtain tempera- 
ture measurements with uncertainties of ~ 15 per cent for X- 
ray luminous galaxy clusters, we used the criterion of ~ 2000 
net counts in the 2-7 keV band per bin (Zhang et al. 2006, 
2007). We therefore require that the screened XMM-Newton 
data of each cluster have more than 4000 net counts in the 2- 
7 keV band after the subtraction of the renormalized blank sky 
background to ensure at least two bins to measure the temper- 
ature profile. Up to the end of 2006, 37 LoCuSS clusters in 
the XMM-Newton archives fulfill this criterion (Fig. I A. 1 1 and 
Table lATT), of which 9 clusters are part of the REFLEX-DXL 
sample (Zhang et al. 2004, 2006; Finoguenov et al. 2005), 12 
clusters are from the LoCuSS Pilot Study sample (Zhang et 
al. 2007), and additionally 6 clusters are from the REXCESS 
sample (Bohringer et al. 2007; Pratt et al. 2007; Croston et 
al. 2007) respectively. Since the REFLEX-DXL, LoCuSS Pilot 
Study, and REXCESS samples are constructed as morphology- 
unbiased samples, our sample, as a collective sample of them 
only with a redshift cut, provides an approximately unbiased 
cluster population. This allows us to investigate (i) the system- 
atics of both X-ray and weak lensing cluster mass estimates, 
and (ii) the morphology dependence of the segregation and 
scatter in the scaling relations. 

All 37 LoCuSS clusters in the XMM-Newton archives were 
uniformly analyzed using the optimized reduction method de- 
veloped for XMM-Newton observations for galaxy clusters at 
redshifts higher than 0.14 (so that the cluster X-ray emission 
covers less than half of the field of view, FOV) with more than 
4000 net photons emitted from the cluster in the 2-7 keV band. 
More details about this method can be found in Zhang et al. 
(2007). The outermost annulus outside the truncation radius 
(8-8.33' or 9.2-10' depending on the truncation radius) was 
used to model the residual background. The XMMSAS v6.5.0 
software was used for data reduction. We present the X-ray de- 

1 http://www.netlib.org/odrpack/ and references therein 



termined cluster center and observational information of the 
sample of all 37 LoCuSS clusters in Table lA.ll 

In our X-ray mass modeling, we used the de-projected 
spectra to measure the spatially resolved radial temperature 
distribution. A double-/? model is adopted to fit the the ICM 
density distribution as follows. Firstly, the slope parameter ft 
is derived by the local fit of the density profile at outer radii 
(> 0.5r5oo) by a single-/? model. As the inner component of the 
ICM density distribution is insensitive to the slope, we fixed 
the slope parameter ft and fit the ICM density distribution to 
derive the inner and outer core radii of a double-/? model. As 
the slope is locally fitted using the outer radii, there should not 
be a significant bias due to the ft steepening in our X-ray mass 
modeling. We therefore only show the mass distribution de- 
rived from the spatially resolved temperature and density dis- 
tributions with no extrapolation to larger radii. With deeper ob- 
servations, the slope of the temperature profile can be better 
determined, which will improve the precision of the mass es- 
timate. The Monte-Carlo error estimation method used for the 
error propagation in our X-ray mass modeling accounts for this. 
More details about the X-ray mass modeling can be found in 
AppendixlAl 

3. X-ray scaling relations: an empirical similarity 

This section describes the detailed modeling and analysis of the 
X-ray data. Readers interested solely in the X-ray /lensing com- 
parison may wish to move straight to the summary in § 13. 131 
and continue from there. 

For this sample the cluster masses are uniformly deter- 
mined using XMM-Newton data. This guarantees minimum 
scatter due to the analysis method. The scaling relations are 
generally parameterized by a power law (Y = FqX 7 ). For each 
relation, we performed the fitting with both normalization and 
slope free. Scatter describes the dispersion between observa- 
tional data and the best-fit values. The scatter in the scaling re- 
lations directly measures the degree of variation of ICM prop- 
erties from the self-similar behavior. We list the best-fit power 
law and the scatter in TableQ] The ODRPACK package weights 
the errors of the data, and estimates the uncertainties of the nor- 
malization parameter. 

To compare the normalization between our sample and the 
samples from other observations and simulations, we also per- 
formed the fitting with the fixed slopes to the values as used for 
their samples. The uncertainty in the difference of the normal- 
ization takes into account the errors of both normalization of 
the two scaling relations for comparison. Using a fixed slope 
might enlarge the segregation when the mass range is not the 
same for two samples. 

Based on the least scattered relations, M^-T and L^-Yx, 
for our sample we checked the self-consistency of the shapes 
(slopes) of the remaining scaling relations. The resulting slopes 
generally appear consistent with the inferred slopes for our 
sample. This allows us to fix the slopes to the self-consistent 
values inferred from our sample and to study the segregation 
and scatter in the scaling relations due to different cluster pop- 
ulations. 
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3.1. Cool cores and X-ray scaling relations 

We investigated the structural similarity of the scaled profiles 
of the X-ray properties in Appendix iBlin order to optimize the 
definition of the cluster global temperature and luminosity to 
obtain least scattered X-ray scaling relations with a brief de- 
scription as follows. The global temperature is derived by the 
volume average of the radial temperature profile limited to the 
radial range of 0.2-0.5r 500 (T .2-o.5r m = T^jrM^ as listed 
in Table IA.2I where Tk and V* are the temperature and vol- 
ume of the A:th radial shell of all the shells (from 1st to nth) 
between 0.2 - 0.5r5oo- We used the values in Table 4 in Zhang 
et al. (2007), and found that the volume average of the radial 
temperature profile gives nearly the same value as the spec- 
troscopic temperature with a spectroscopic temperature versus 
volume averaged temperature ratio of 0.98. According to sim- 
ulations from Rasia et al. (2005), our temperature definition 
will be lower than the emission measured temperature by 20- 
30 per cent. The core radii populate a broad range of values 
up to 0.2 rsoo. With the bolometric X-ray luminosity includ- 
ing and excluding the < 0. 2^500 region (Table IA.2I ). we show 
the normalized cumulative cluster number count as a function 
of the fraction of the luminosity attributed to the cluster core 
(< 0.2r5oo) in Fig. IB. 71 Up to 80 per cent of the bolometric 
X-ray luminosity is contributed by the cluster core (< 0.2r5oo). 
We therefore use the X-ray luminosity corrected for the cool 
core for the R < 0.2r5oo region (L^ in Table lA.2l . 

3.2. Mass versus temperature 

Hereafter Mf Q[j denotes the total cluster mass derived from X- 
ray data at a density contrast of 500. In Fig. [2] we present the 
M* Q() -T relation for our sample and its best fit, which gives a 
slope of 1.65 +0.26. 

We also performed the fitting with the slope fixed to the 
values 1 .49 and 1 .489 (consistent with our fit 1 .65 + 0.26) as 
used in Arnaud et al. (2005) and Vikhlinin et al (2006a). Our 
M^ Q() -r relation agrees (1+7 per cent) with the M-T rela- 
tions for 13 low-redshift clusters (z < 0.23, 0.7-9 keV) using 
Chandra observations (Vikhlinin et al 2006a), and for 6 relaxed 
nearby clusters (z < 0.15, T > 3.5 keV) using XMM-Newton 
data (Arnaud et al. 2005), respectively. No evident evolution is 
found for the M-T relation after the redshift evolution correc- 
tion due to LSS growth (defined in AppendixlBl. 

Under the assumption that the temperature is constant 
within the cluster, the Mf 0Q -T relation based on ROSAT ob- 
servations and the temperature measurements from literature 
related to the HIFLUGCS sample (Reiprich & Bohringer 2002; 
Popesso et al. 2005; Chen et al. 2007) gives higher normaliza- 
tion by ~ (40-50) + 20 per cent than the Mf^-T relation based 
on Chandra/XMM-Newton observations. As the temperature 
distribution is used in the mass modeling to obtain the XMM- 
Newton/ Chandra results, the high normalization in the ROSAT 
results is most probably due to the isothermal assumption in 
the mass modeling. We fixed the slope to 1.48 as used for the 
hot clusters (> 3 keV) with ASCA measured temperature pro- 
files in Finoguenov et al. (2001), and found the normalization 
is higher by 19 + 21 per cent than the value in Finoguenov et 



al. (2001). It indicates the importance of well defined global 
temperature based on a spatially resolved temperature profile 
on the scaling relation studies. 

We also fixed the slope to 1.59 (consistent with our fit 
1 .65 + 0.26), and found that the normalization of our X-ray M- 
T relation is lower by 19 + 6 per cent than the one from simu- 
lations (SPH code) including radiative cooling, star formation, 
and SN feedback in Borgani et al. (2004). This is similar to the 
finding in Nagai et al. (2007a) that the X-ray measured clus- 
ter masses within rsoo are lower than the cluster masses from 
simulations. 

The M-T relation of our observations uses the temper- 
ature comparable to the spectral-like temperature. The M-T 
relation from simulations in Borgani et al. (2004) uses the 
emission-weighted temperature. As Rasia et al. (2005) men- 
tioned that the normalization of the M-T relation from sim- 
ulations in Borgani et al. (2004) is higher by ~ 30 per cent 
using the spectral-like temperature instead of the emission- 
weighted temperature, 19 + 6 per cent gives the minimum dif- 
ference between our observations and simulations in Borgani 
et al. (2004). 

Evrard et al. (1996, emission weighted temperature) sim- 
ulated ROSAT observations of 58 nearby clusters (z ~ 0.04, 
1-10 keV), for which the normalization of the M-T relation 
is higher by 40 ± 16 per cent than ours (with the slope fixed 
to 1.5). It might be due to the neglect of complex gas physics 
(e.g. radiative cooling, star formation and SN feedback) in their 
simulations. 

3.3. Mass versus Y x 

Simulations (Kravtsov et al. 2006) indicate that the X-ray ana- 
log of the integrated SZ flux (Fx — ^gas • T) can be used as a 
robust low scatter cluster mass indicator with less scatter com- 
pared to other X-ray observables. Kravtsov et al. (2006, Fig. 7 
in their paper; also see Table 4 in Nagai et al. 2007b) found 
the M-Yx relation for the sample in Vikhlinin et al. (2006a) 
confirms the low scatter and shows an empirical similarity. 
Maughan et al. (2007) investigated the M-Yx relation using 
Chandra observations and found that the scatter varies with 
the data quality. It is therefore important to calibrate the M-Yx 
relation and to check whether Yx does show the least scatter 
as a mass proxy using high quality observations from XMM- 
Newton. 

We present the M-Yx relation for the sample in Fig. [T] 
in which the best fit gives a slope of 0.624 ± 0.061, consis- 
tent with 0.526 + 0.038 derived from Chandra data (see Fig. 7 
in Kravtsov et al. 2006 and Table 4 in Nagai et al. 2007b). 
We fixed the slope to 0.526, and found an agreement (2+12 
per cent) of the normalization between M-Yx relations for our 
sample and theirs. With the slope fixed to 0.548, the normaliza- 
tion for our sample is higher than for the sample in Arnaud et 
al. (2007) by 6 per cent, which is also based on XMM-Newton 
observations. There is no obvious evolution for the M-Yx rela- 
tion after the redshift evolution correction due to LSS growth. 
We found that the scatter in Yx for M-Yx is similar to the scat- 
ter in T for M-T. The comparison between our sample and the 
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samples in Vikhlinin et al. (2006a), Maughan et al. (2007) and 
Arnaud et al. (2007) might indicate that the slope of the M-Yx 
relation varies as a function of the mass range of the sample. 

As Y x appears to be a good mass proxy, the cluster mass can 
be estimated from Y x via the M-Y x relation. Hereafter such a 
mass estimate is called the Y x inferred cluster mass, M Yx . 

We also fixed the slope to 0.568 (consistent with our fit 
0.624 + 0.061) as obtained in simulations (Nagai et al. 2007b), 
and found the normalization of our X-ray M-Y x relation is 
lower than the one in simulations from Nagai et al. (2007b, the 
mass-observable relation using the spectroscopic temperature 
and true mass; also see Kravtsov et al. 2006) including cool- 
ing and galaxy feedback by 18+4 per cent. From simulations 
in Nagai et al. (2007b), the hydrostatic mass estimate is about 
20% lower than the true mass, similar to the differences we 
observe. Here we compare our mass measurements with "the 
truth" from simulations rather than the hydrostatic masses de- 
rived from the simulations in Nagai et al. (2007b). In part this 
is because it gives a consistent basis upon which to compare 
lensing and X-ray masses with the simulations. 

3.4. Mass versus gas mass 

The total mass versus gas mass relation (M-M gas ) for our sam- 
ple is shown in Fig.Q] giving a slope of 0.906+0.078, consistent 
with 0.81 1 + 0.067 derived from Chandra data (Kravtsov et al. 
2006, Nagai et al. 2007b). We fixed the slope to 0.811 as de- 
rived for the sample in Kravtsov et al. (2006), and found that 
the normalization of the M-M gas relation for our sample agrees 
(2+11 per cent) with theirs. No evident evolution is found for 
the M-Mg as relation after the redshift evolution correction due 
to LSS growth. 

We also fixed the slope to 0.894 (consistent with our fit 
0.906 ± 0.078) as derived in Nagai et al. (2007b), and found the 
normalization of our X-ray M-M gas relation is lower than the 
one in simulations in Nagai et al. (2007b) by 13+4 per cent. 

3.5. Gas mass versus temperature 

Fig.[3]shows the M gas ^oo-T relation for our sample and its best 
fit giving a slope of 1.859 ± 0.187. 

We obtained a slope of 1.84 for the M ga ^0{)-T relation for 
the sample in Vikhlinin et al. (2006a) using the scaling relations 
in Kravtsov et al. (2006). We then fixed the slope to 1.84 (con- 
sistent with our fit 1.859 + 0.187), and found that the M gaSj50 o- 
T relation for our sample agrees (3 + 11 per cent) with their 
nearby relaxed cluster sample based on Chandra observations. 
No evident evolution is found for the M gas -T relation after the 
redshift evolution correction due to LSS growth. 

Castillo-Morales & Schindler (2003) investigated a sample 
of 10 nearby clusters (0.03 < z < 0.09, 4.7-9.4 keV) observed 
by ROSAT and ASCA. The normalization for their samples is 
higher by 19 + 3 per cent (no errors available for the normaliza- 
tion in Castillo-Morales & Schindler 2003) than ours with the 
slope fixed to 1.8 as used in their work. The M gas ,500-r relation 
in Mohr et al. (1999) is based on a flux-limited sample of 45 
nearby clusters spanning a temperature range of 2-10 keV also 



observed by ROSAT and ASCA. Their normalization is higher 
than ours by 81 ±7 per cent using their slope of 1.98. Chen 
et al. (2007) present the M gaS 5oo-r relation based on the non- 
cool core clusters (non-CCCs) with T > 3 keV in HIFLUGCS, 
in which the normalization is higher than ours by a factor of 
2 with the slope fixed to 1.8. The difference in the normal- 
ization might be due to the difficulty of performing spatially 
resolved temperature measurements with ROSAT/ASCA data. 
At that time, one had to use the rsoo defined from simulations 
(e.g. Evrard et al. 1996) which is larger than the value obtained 
from observations based on the radial temperature profile using 
XMM-Newton data. 

We also fixed the slope to 1.8 (consistent with our fit 
1.859 + 0.187) as shown in Borgani et al. (2004), and found 
the normalization of our X-ray M g!ls -T relation is lower than 
the one from simulations in Borgani et al. (2004, emission 
weighted temperature) by 7 + 1 1 per cent. If only consider- 
ing the massive systems, we observe a slightly smaller segre- 
gation between our observations and their simulations as shown 
in Fig. [3] 

3.6. Luminosity versus Y x 

The detailed calculation of the X-ray luminosity is described in 
Appendi) jB.4l The Lboi-lx relation for this sample is shown in 
Fig. [3] with a best-fit slope of 0.951 + 0.081. 

We fixed the slope to 0.94 (consistent with our fit 0.95 1 + 
0.081) as used in Maughan (2007), and found that the normal- 
ization for our sample is higher than for the sample of 1 16 mas- 
sive clusters (0.1 < z < 1.3) in their work by 10 ± 4 per cent. 
Compared to the observational errors, there is no obvious evo- 
lution in the Lboi-lx relation after the redshift evolution cor- 
rection due to LSS growth. 

3.7. Luminosity versus temperature 

The L-T relations and their best fits, giving slopes of 2.127 + 

0.323 for L c °^ 24 -T and 2 - 614 ± 0319 for L b°"- r ' are shown 
in Fig. g] 

The L-T relation has been thoroughly investigated for 
nearby cluster samples based on ROSAT and ASCA obser- 
vations, for example, in (i) Markevitch (1998, 30 clusters, 
0.04 < z < 0.09, T > 3.5 keV), (ii) Arnaud & Evrard (1999, 24 
clusters, z < 0.37, T > 2 keV), (iii) Reiprich & Bohringer 
(2002, in which the flux-limited morphology-unbiased sam- 
ple HIFLUGCS was initially constructed and investigated. The 
fits used here are for the whole HIFLUGCS sample includ- 
ing groups), (iv) Ikebe et al. (2002, a flux-limited sample of 
62 clusters, z < 0.16, 1-10 keV), and (v) Chen et al. (2007, 
HIFLUGCS). 

The best-fit slope 2.614 + 0.319 for the L™"-T relation for 
our sample rules out the standard self-similar prediction L ozT 2 
at 2cr significance. The normalization for our sample agrees 
with the normalization for the nearby samples using represen- 
tative non-CCCs in Markevitch (1998) on 6 + 1 1 per cent level 
with the slope fixed to 2.64, in Arnaud & Evrard (1999) on 
4 + 9 per cent level with the slope fixed to 2.88, and in Chen et 
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al. (2007) on 13+22 per cent level with the slope fixed to 2.78, 
respectively. 

The normalization of the L^™_ 2 4keV -T relation for our sam- 
ple is slightly lower than for the nearby samples in Markevitch 
(1998) by 13 ± 11 per cent with the slope fixed to 2.10, but 
shows an agreement with the sample in Ikebe et al. (2002) on 
5 ± 38 per cent level with the slope fixed to 2.44. 

The observed evolution in the L-T relation is consistent 
with the prediction of the LSS growth and is therefore not seen 
on the corrected relations. 

There is no significant difference between the normaliza- 
tion of the L-T relations using ROS AT/ASCA data and XMM- 
Newton/ 'Chandra data, which is not the case for the M-T and 
M g . ds -T relations. The spatial resolution is of prime importance 
to get precise temperature gradient for a good X-ray mass mea- 
surements. However, L and T can be derived reliably without 
spatially resolved spectroscopy. Theforefore the L-T relation 
from ASCA/ROSAT is more likely to be in agreement with 
Chandra/XMM than mass-observable relations. 

The normalization is higher by 41 ±24 per cent using a fixed 
slope of 2.598 for Lg°£_ 2 4teV -T and by 36 ± 23 per cent with 
the slope fixed to 2.98 for Lg^-T for the HIFLUGCS sample 
in Reiprich & Bohringer (2002). The temperature can be biased 
low and the X-ray luminosity can be biased high for CCCs in 
the HIFLUGCS sample. At the same time both temperature and 
luminosity can be boosted high for merging systems. It is dif- 
ficult to predict the combined effects without looking into the 
detailed temperature and density profiles for the HIFLUGCS 
sample. 

3.8. Luminosity versus mass 

The L£°"-M relation is shown in Fig. [2] giving a best-fit slope 
of 2.325 +0.701. The scatter in the L-M relation is much larger 
than the scatter in the other X-ray scaling relations for our sam- 
ple. 

The normalization of the L-M relations for our sam- 
ple is higher than for the HIFLUGCS sample based on 
ROSAT/ASCA observations (Reiprich & Bohringer 2002; 
Popesso et al. 2005; Chen et al. 2007) by a factor of ~ 2. The 
difference in the normalization might be due to the isothermal 
assumption in their mass modeling and boosting effects due to 
CCCs in HIFLUGCS. 

We also fixed the slope to 1.572 (smaller than our fit 
2.325 ± 0.701) as derived from the Lb \-T and M-T relations 
in Borgani et al. (2004, emission weighted temperature), and 
found the normalization of our X-ray Lb \-M relation is lower 
than the one from simulations in Borgani et al. (2004) by 18+32 
per cent. 

3.9. r 50 o versus global temperature 

The best fit for our sample gives r^m = 
10 -o.38±o.07 r o.60±o.08^ (z:) -i Mpc similar t0 the finding in 

Arnaud et al. (2005) observing r 5Q0 oc r°-5°±o.05 for 6 re _ 
laxed nearby clusters (z < 0.15, 2-9 keV, T > 3.5 keV) 
using XMM-Newton data. Pacaud et al. (2007) found 



r 50 o = 0.375(77/fceV) - 63 /j7 3 (z)- 1 Mpc based on the measure- 
ments in Finoguenov et al. (2001), which is in agreement with 
ours within the uncertainty. 

3.10. Self-consistent scaling relations 

As shown in Table [T] both M gas -T and Lboi-lx relations show 
smaller scatter than the remaining scaling relations. The er- 
rors (10 per cent) of the slope parameters are also the small- 
est for these 2 relations compared to the other relations. We 
therefore checked the self-consistency of the empirical scaling 
relations based on these 2 relations (M gas oc t , 1-859±o.187 an( j 
Lbo\ K y°- 951±0 - 081 ) for our sample as follows. 

The definition of the X-ray bolometric luminosity gives 
Ami K M^r^Af -1 . Combining with L boX oc T°- 951 , we derive 
M 2 m T°' 5 M- 1 = (M gas • r) 951 . Given M gas oc T 1859 , the M-T 
relation is M oc T 1A ", and the M-M gas relation isMoc M° a 806 , 
respectively. 

Combining M gas oc r 1859 and M oc T lA " yields M oc 

y0.524 
I X 

We obtain Lboi K T 2 J19 on the basis of L b oi K T° 951 and 
M gas oc r L859 , and L bol oc M 1814 on the basis of L bol oc r 2 ' 719 
and M oc T 1A ", respectively. 

The slopes of the empirical scaling relations for our sample 
inferred from M gas oc r 1 - 859 * - 256 and L bol oc y°- 9 5i±o.08i are in 
agreement with the best-fit slope values in Table Q] This indi- 
cates that our sample shows a self-consistent, empirical simi- 
larity in the scaling relations. 

3.11. Cool core dependence of the normalization and 
scatter 

Different cluster populations might show different normaliza- 
tion and scatter in the scaling relations. Most CCCs show sim- 
pler morphologies than non-CCCs, though this is not always 
the case, e.g. Abell 2390 (Allen et al., 2001; Swinbank et al., 
2006). We checked the cool core dependence of the normal- 
ization and scatter using our sample. There is no significant 
bimodal distributions in the X-ray scaling relations. However, 
the CCCs and non-CCCs defined in Appendix lAl are still used 
to divide the sample into two subsets in order to quantify the 
segregation in the scaling relations due to cool cores. We fixed 
the slope parameters to the prediction in § 13.101 on the basis 

Of M m oc 7-l-859±Q.256 an( j ^ x ^0.951*0.08^ an( j compare d 

the normalization and scatter in the scaling relations between 
the whole sample and non-CCCfl (Table Q]). The normaliza- 
tion agrees to better than 17 per cent. The difference in scatter 
is most marked in the L-M relations (also cf. Chen et al. 2007). 
This may be due to the use of the hydrostatic equilibrium as- 
sumption for the determination of the X-ray cluster mass. 



2 The CCC subsample only consists of 9 clusters as shown in 
Table IbTI which is statistically less sufficient compared to the non- 
CCC subsample to investigate the scatter, specially for the mean de- 
termined using a Gaussian model fit to the scatter histogram. 
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3. 12. Mass-observable relations from X-ray 
observations and simulations 

In the section we summarized the comparison between our 
sample and the samples from numerical simulations. Given a 
wide variety for sources of the simulated data, we provide here 
a comparison between the different methods and also comment 
on the sensitivity of various predictions in the adopted schemes. 
In simulations the measurements, such as the total mass and gas 
mass, are known precisely, while some measurements, such as 
the X-ray luminosity, are calculated in a straightforward man- 
ner, and some measurements, such as the observed tempera- 
ture, are method-dependent. For example, Rasia et al. (2005) 
mention that the normalization of the M-T relation from simu- 
lations in Borgani et al. (2004) is higher by ~ 30 per cent using 
the spectral-like temperature instead of the emission weighted 
temperature. We used the following simulations in the com- 
parison: Evrard et al. (1996, emission weighted temperature), 
Borgani et al. (2004, emission weighted temperature), and 
Nagai et al. (2007b, spectroscopic temperature). Further details 
regarding the numerical resolution, cosmological parameters, 
and codes (grid versus particle) of simulations can be found in 
the papers cited above. 

The normalization of the X-ray mass-observable relations 
is lower than the one from simulations (M-Yx and M-M gas in 
Nagai et al. 2007b; M-T and L hol -M in Borgani et al. 2004, 
emission weighted temperature) by(l) 18+4 per cent for M- 
Y x , (2) 13 +4 per cent for M-M gas , (3) 19 + 6 per cent for M-T, 
and (4) 18 + 32 per cent for L\, \-M, respectively. 

The upward correction for the X-ray masses was suggested 
to be due to the steepening of the ICM density and tempera- 
ture distributions assuming hydrostatic equilibrium by the com- 
parison between Chandra observations (Kravtsov et al. 2006) 
and simulations (Nagai et al. 2007b). For our sample, the den- 
sity profiles are well resolved spatially (> 100 bins) while the 
temperature profiles are poorly resolved spatially (> 2 bins). 
Luckily, the universal temperature profile in Vikhlinin et al. 
(2005) is based on the temperature profiles well resolved spa- 
tially for 13 low redshift galaxy clusters. We therefore investi- 
gated the steepening of the ICM density and temperature pro- 
files in the outer region using our sample for the density distri- 
bution, and the sample in Vikhlinin et al. (2005) for the tem- 
perature distribution. We found that the mass can be higher by 
up to 25 per cent at r^m taking into account our ICM density 
steepening and the temperature steepening using the universal 
temperature profile in Vikhlinin et al. (2005) in the hydrostatic 
equilibrium assumption. One should be cautious with this point 
since we have almost no resolution on temperature profiles be- 
tween 0.5r5oo and rsoo, where the temperature profile steepen- 
ing can not be verified. One requires deep X-ray observations 
with 8-10 radial bins for the temperature measurements on the 
per cent precision level covering the area of more than 0.75r5oo 
for a representative sample of more than 30 clusters to suffi- 
ciently verify the steepening of the ICM temperature profiles 
to investigate the bias. 



3.13. Summary of the X-ray scaling relations 

The above comparison of the X-ray scaling relations between 
this sample and published nearby and more distant samples 
shows that the evolution of the scaling relations can be ac- 
counted for by the redshift evolution due to LSS growth as 
given in Appendix iBl 

Based on our sample, the scatter in the X-ray mass- 
observable scaling relations is (0.03-0. 15dex). 

The normalization agrees to better than 17 per cent for the 
scaling relations between the sample and the non-CCC sub- 
sample, with by far the worst case being the L^-M relation. 
In general, the slopes of the scaling relations (e.g. 1.859+0.187 
for M gas -r) indicate the need for non-gravitational processes. 
This fits a general opinion that galaxy clusters show an empiri- 
cal similarity. 

The normalization of the X-ray mass-observable relations 
(i.e. M-T, M-Yx and M-M gas ) for our sample is lower by 13- 
19 per cent with less than 6 per cent uncertainties than the re- 
lations from simulations (Borgani et al. 2004; Kravtsov et al. 
2006; Nagai et al. 2007b). 

4. rj *' x and the reduced scatter in the 

mass-observable relations and L bo] -Y x relation 

As shown in the left panels in Figs.[TJ{2l a few merging systems 
show significant deviation from the best-fit mass-observable 
relations. Since simulations (Nagai et al. 2007b) have shown 
that Yx is claimed to be a good mass proxy even for merg- 
ing clusters in which hydrodynamical equilibrium may be dis- 
turbed, we may improve the radial scaling by using a mass pro- 
file based on Y x as a mass estimator. We used the following 
method to reduce the deviation in mass-observable relations 
for merging clusters by using the measured masses and X-ray 
observables at a new radial scale, r Y ^ x . 

The radial scale r™ is derived as follows. The procedure 
is initiated by deriving the best-fit M-Yx relation using the X- 
ray masses M500 and Yx parameters at the X-ray determined 
^500 for the 37 LoCuSS clusters with the slope fixed to 0.568, 
the prediction from simulations in Nagai et al. (2007). This best 
fit is used to calculate the mass profile M Yx (r) from the profile 
Yx(r) = M gas (r) * 7o.2-o.5r 50 o ■ The radial scale, within which the 
density contrast derived from the mass profile M Yx {r) is 500, 
is called r test . 

We re-fit the M-Yx relation but using the X-ray masses and 
Yx parameters at r test for the 37 LoCuSS clusters with the fixed 
slope of 0.568. This best fit is used to re-derive the mass profile 
M Yx (r) from the profile Yx(f), and M Yx (r) can then be used 
to re-determine the radial scale r test , within which the density 
contrast, derived from the mass profile M Yx (r), is 500. We iter- 
atively re-fit the M-Yx relation using the X-ray masses and Yx 
parameters at r test , and re-determine the radial scale r test using 
M Yx (r) till the best fit of the M-Yx relation does not vary any 
more for our sample. The iteration is less than 10 times to ful- 
fill the goal. The radial scale, J*™ , is the final r test . The radial 
scale r Y y' x and the X-ray mass and Yx parameter at r™ for 
each cluster are listed in Table IC.ll 
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In Figs. [5}|6] we show the mass-observable scaling rela- 
tions using the X-ray masses and observables at r 5( 5, for our 
sample. The deviation from the best fit is indeed reduced for 
merging systems. 

The normalization of the mass-observable relations using 
X-ray masses and observables measured at rj *' x is lower than 
the one from simulations by, (1) 24 ± 3 per cent for M-Yx, (2) 
17 +4 per cent for M-M gas , (3) 17 + 5 per cent for M-T, and (4) 
17 + 31 per cent for L^-M, respectively. The difference in the 
normalization from X-ray observations and simulations does 
not vary significantly compared to the investigation using the 
X-ray masses and observables measured at rsoo- This difference 
in the normalization is relatively independent of the scatter in 
the X-ray scaling relations. 

5. Weak lensing versus X-ray mass comparison 

Weak lensing and X-ray approaches provide independent clus- 
ter mass estimates. Comparison of X-ray and weak lensing 
masses is a powerful test of systematic errors in cluster mass 
determination. In total, weak lensing mass estimates are avail- 
able in the literature for 19 of the total sample of 37 clus- 
ters. Ten of the 19 were observed over a wide-field (i.e. out to 
R > r 20 o) with CFHT by Bardeau et al., (2005; 2007 - hereafter 
the B07 sample). 

Fifteen of the 19 were studied by Dahle (2006 - here- 
after the D06 subsample) using the shear estimator based on 
the method in Kaiser (2000) but applied for small-field (2k 
chip, typically to 0.4r 50 o) data from NOT and UH88. For this 
shear estimator, the polarizability includes calculating vari- 
ous derivatives with respect to galaxy flux, size and elliptic- 
ity. In practice, this is done numerically by binning the ob- 
served galaxies in a 3-D flux-size-shape space. For small-field 
data from NOT and UH88, there are only ~ 10 3 suitable back- 
ground galaxies available for each cluster, which gives sparsely 
populated bins and thus introduces some additional uncertain- 
ties to the calculation of the polarizability. This shear estimator 
overcomes the problem by including galaxy information from 
a number of different fields before the binning and polarizabil- 
ity calculation, and produces more stable results. The overlap 
between the D06 and B07 samples is 6 clusters, giving a net 
literature-based weak lensing sample of 19 clusters. 

The following sections describe the mass comparison using 
M500 (§ 15. It , and the mass measured at rj *' x (§ I5.21 i. 

5.1. M500 comparison 

To avoid the uncertainties caused by the extrapolation of the 
X-ray measured mass profile beyond rsoo in the projection, 
we present the weak lensing versus X-ray mass calibration us- 
ing spherical masses within radii of a density contrast of 500 
(M500). Hereafter M^ denotes the weak lensing mass within 
weak lensing measured rsoo. The average ratio of M™ ' () to Mf QQ 
gives 1.7 + 0.2 for the D06 subsample, and 1.1 + 0.2 for the 
B07 subsample, respectively. The better agreement obtained 
from B07 subsample arises from the large field of view of the 
CFHT data used by B07 (see above), in contrast the D06 data 



do not reach rgoo, and thus the D06 masses require extrapola- 
tion. This shows the advantage of the LoCuSS project using 
uniformly observed and reduced data with large FOV to cali- 
brate the weak lensing and X-ray masses and to investigate the 
scaling relations and their scatter. 

The parameter Yx shows less scatter than the X-ray mass 
as shown in the M-Yx relation. To check whether the scatter in 
the weak lensing versus X-ray mass comparison is due to the 
scatter in X-ray masses, we replace the X-ray mass estimates 
by the mass given by Yx at rgoo and the X-ray Mgoo-Jx relation, 
Mj* = M(Y x (r 500 )). 

In the upper panels in Fig. [7] we show M^J against M^, 
and M 5Q q, respectively, for the D06 and B07 subsamples. The 
scatter is reduced using compared to the scatter using 
Mf m . The slope of the best-fit power law is 1.61 + 0.39 Of 2 = 
0.06/23) for the M™' versus Mf m relation, and 1.57 + 0.38 
(X 1 = 21.2/23) for the M™' () versus M^ X relation. The reduced 
X 2 for the M^ Q versus Mf m relation is much smaller than 1, 
which might indicate the over-estimation of the errors of Mf Q(j . 
When the slope is fixed to 1, the normalization of the best-fit 
power law is 1.28 + 0.14 for the M^ versus relation, 
and 1.21+0.13 for the MfL versus mJ Q q relation, respectively. 
On average, the weak lensing masses (M^L) are higher than 
the X-ray masses (Mf QQ ) by 28 + 14 per cent according to the 
normalization. Systematically higher weak lensing masses than 
X-ray masses at rgoo have also been reported by Mahdavi et al. 
(2007). 

The weak lensing to X-ray mass ratio derived above takes 
into account the error weighting in both X-ray and weak lens- 
ing masses. The error in X-ray masses does not vary much from 
cluster to cluster. However, for a few clusters, the error in the 
weak lensing mass is a few times as large as for the remaining 
clusters. This may introduce a bias in the derived weak lensing 
to X-ray mass ratio using the above error weighted method. 
To avoid such a bias, we checked the mean of the M^ to 
M^ )() ratios and M™ ' to M 50 g ratios as follows. In the lower 
panels in Fig. [7] we show the normalized differential cluster 
number count as a function of the M500 ratio. It can be fit- 
ted by a Gaussian model Y = Y exp- ((x - Xo),w)2/2 . The his- 
togram of M^/M^ gives F = 0.15 +0.01, X = 0.99 + 0.05, 
and W = 0.51 ± 0.05. And the histogram of M^ /M^ gives 
To = 0.15 + 0.01, X = 1.00 + 0.05, and W = 0.45 ± 0.05. 
In logarithmic space, the best Gaussian fit provides mean val- 
ues lower than the values obtained in linear space by ~ 15 per 
cent for both cases. In linear space, the mean is 0.99 with 51 
per cent scatter in M^/Mf m , and 1.00 with 45 per cent scat- 
ter in M^/M^, respectively. As shown in the lower part of 
the lower panels in Fig. [7] the discrepancy between and 
MS> (or M^ q) can be up to 5 times the errors of M%L for in- 
dividual clusters. To our surprise, the clusters in the tails of 
the Gaussian distribution show both pronounced mass discrep- 
ancies and very small errors in weak lensing masses. It indi- 
cates that the clusters showing large mass discrepancies may 
be over-weighted in the determination of the weak lensing to 
X-ray mass ratio using the error weighted method. As a conse- 
quence, the deviation of the weak lensing to X-ray mass ratio 
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from 1 can be biased high using the error weighted method. 
On the contrary, the Gaussian fit to the histogram of the weak 
lensing to X-ray mass ratio gives the value of the dominant 
population. 

In the above test, we used the best-fit X-ray M-Yx relation 
for our sample with the slope fixed to 0.568 predicted from 
simulations (Nagai et al. 2007b) to derive Mj x from Ix.soo- 
We made an additional attempt to repeat such tests by using 
the best-fit X-ray M—Yx relations for our sample with (i) a free 
slope, and (ii) the slope fixed to 0.524 as predicted from the 
empirical self-similarity of our sample described in § 13.101 re- 
spectively. The mean of the MfL to ratios varies within 
14 per cent, and the scatter varies within 7 per cent in the tests 
via different X-ray M-Yx relation. Within the scatter, it holds 
that the mean of the weak lensing to X-ray M500 ratio is about 
1. 



5.2. Mass comparison at r^ x ' x 



500 

In the above mass comparison, we used M500 for X-rays taken 
at the X-ray determined r 5 oo and M500 for weak lensing taken 
from the weak lensing determined r 5 oo. This can introduce two 
sources of difference, (1) the normalization of the mass profile, 
and (2) the radius up to which the mass is integrated. To reduce 
such differences, we investigated the mass comparison using 
the masses measured at the same radius as follows. 

Despite the mass comparison using both weak lensing 
masses and X-ray masses measured at the radial scale of r Y x ' x , 
we also made a comparison between the weak lensing mass es- 
timates and the mass given by Yx at r™ X and the X-ray M500- 



Y x relation, M Fx (r^ x ' x ) 



M(Y x (r^ x )), as shown in Fig. [8] 



'500 ' "" V " AV '500 

When the slope is fixed to 1, the normalization of the best- 
fit power law for the M wl (r^ x ' X ) versus M x (r^ x ) relation is 
1 .09+0.08, and for the M wl (rJ x ; x ) versus M Yx (rj* ,x ) relation is 
1 .04 + 0.08, respectively. On average, the weak lensing masses, 
M wl (rJ ( x ' x ), are higher than the X-ray masses, M x (rJ x ' x ), by 
9 + 8 per cent according to the normalization. Following the ar- 
guments presented in Churazov et al. (2007), the current com- 
parison limits the non-thermal pressure support to 9 ± 8 per 
cent according to the ratio between the weak lensing and X-ray 
mass estimates. 

With a Gaussian fit Y = Y exp- ({x - x o )/W)2/2 , the his- 
togram of M wl (rJ x ' x )/M x (rJ x ' x ) gives Y Q = 0.19 ± 0.02, 
X = 1.16 ± 0.04, and W = 0.42 + 0.04. And the his- 
togram of M wl (rJ x ' x )/M Yx (rJ () x ' x ) gives Y = 0.22 ± 0.02, 
X = 1.01 ± 0.03, and W = 0.33 + 0.03. The scatter is there- 
fore reduced to 33^+2 per cent compared to 45-51 per cent for 



almost independent of X-ray mass, (3) global temperature, 
and (4) X-ray bolometric luminosity. We present the mass- 
observable relations using weak lensing masses for the subsam- 
ple of the 19 clusters on the right column of panels in Figs.[T]-[2] 

6. 1 . Scatter and systematics 

The scatter using weak lensing masses is larger than the scatter 
using X-ray masses by about a factor of 2 (Table [2]). This dif- 
ference is probably caused by a combination of (i) larger mea- 
surement errors on the weak lensing mass measurements than 
on the X-ray mass estimates, and (ii) possible under-estimation 
of the scatter in the X-ray scaling relation due to intrinsic cor- 
relation between the axes when both are based on the same 
dataset. 

We also seek a physical origin of the scatter in the lens- 
ing relation, by dividing the clusters into single and merger X- 
ray morphologies (defined in Appendix lA.2b . The single clus- 
ters only show significantly less scatter in mass at fixed ob- 
servable than the merger clusters in the luminosity-mass rela- 
tion, the reduction in scatter being ~50%. This indicates that, 
in general, merging systems introduce insignificant scatter in 
the weak lensing mass based mass-observable relations at rsoo- 
The larger scatter in the luminosity-mass relation of merger 
clusters may indicate shock heating of gas in merger clusters, 
although it is odd that this effect does not also affect the mass- 
temperature relation. There is also no evidence for a difference 
in normalization between single and merger clusters in any of 
the scaling relations examined above except for the L\, \-M re- 
lation, showing 32 per cent difference. This might indicate that 
boost effects are significant on the X-ray luminosity than on the 
temperature for merging systems. 

The scatter in the scaling relations is nearly symmetric us- 
ing weak lensing masses irrespective of being classified as 
CCCs or non-CCCs (Figs. QJO. The normalization is nearly 
the same using weak lensing masses between all 19 clusters and 
the non-CCC subsample, giving up to 7 per cent difference (i.e. 
7 per cent difference in the L\, \-M relation). This indicates low 
systematics in the weak lensing mass based mass-observable 
scaling relations concerning being CCCs and non-CCCs. 

6.2. Observations versus simulations 

To study the difference in the normalization between the mass- 
observable relations using weak lensing masses and the one 
from simulations, we fixed the slope to the value derived from 
simulations and fitted the normalization for the 19 LoCuSS 
clusters. 



6. Weak lensing mass M™, based 
mass-observable relations 

We now investigate the scaling relations between weak lens- 
ing mass and four representative X-ray observable proxies of 
the cluster mass: (1) Yx, which has been suggested as the 
mass proxy with least scatter on the basis of simulations (e.g. 
Kravtsov et al. 2006), (2) gas mass, which is an observable 



6.2.1. First Attempt 

First we use weak lensing masses computed within rsoo derived 
from the lensing data with X-ray observables measured at the 
X-ray determined rsoo- 

We fixed the slope to 0.568 as obtained from simulations 
(Nagai et al. 2007b) for the M-Yx relation, and found that the 
normalization using weak lensing masses is marginally lower 
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than the one from simulations (AMR code) in Nagai et al. 
(2007b, also see Rravtsov et al. 2006) by 1 + 10 per cent. 

With the slope fixed to 0.894 predicted from simulations 
(Nagai et al. 2007b) for the M-M gas relation, we found a 
marginally higher normalization of the M-M gas relation us- 
ing weak lensing masses than the one from simulations (AMR 
code) in Nagai et al. (2007b) by 3 + 11 per cent. 

Using the slope fixed to 1 .59, we found that the normaliza- 
tion of the weak lensing mass based M-T relation is marginally 
higher by 4 + 12 per cent than the one from simulations (SPH 
code) including radiative cooling, star formation, and SN feed- 
back in Borgani et al. (2004). 

The Lb \-M relation from simulations can be derived from 
the M-T and Lboi-r relations from simulations (SPH code) in 
Borgani et al. (2004), which shows a slope of 1.572. With the 
slope fixed to 1.572, we found the normalization of the L\, \-M 
relation using weak lensing masses is lower than the one from 
simulations by 52 ± 34. 

In summary, there is no statistically significant difference 
between the normalization of the observed and simulated rela- 
tions when using lensing mass and X-ray observables derived 
from within r^m obtained from each respective dataset. 



6.2.2. All measurements within rj *' x 



In § 16.2.11 the scatter in the weak lensing mass based scal- 
ing is so much larger than the difference in the normalization 
between observations and simulation. It is therefore not clear 
whether there is consistency between weak lensing mass based 
mass-observable relations and mass-observable relations from 
simulations. In Table [3] we present the best-fit normalization 
of weak lensing mass based mass-observable relations using 
the measurements at r™ listed in Table IC.ll As shown in 
Figs. l9TfT0l the scatter is greatly reduced compared to the re- 
lations using the weak lensing masses M" ' () at the weak lensing 
determined rsoo and X-ray observables measured at the X-ray 
determined 7-500. 

As derived in §|4] the normalization of the mass-observable 
relations using X-ray masses and observables measured at 
r 50o' X is lower than the one from simulations by, (1) 24 ± 3 per 
cent for M-Yx, (2) 17 + 4 per cent for M-M gas , (3) 17 + 5 per 
cent for M-T, and (4) 17 + 31 per cent for L\, \-M, respectively. 
The normalization of the mass-observable relations using the 
weak lensing masses and X-ray observables measured at r^™ 
is lower than the one from simulations by, i.e. (1) 18 + 8 per 
cent for M-Y x , (2) 17 + 8 per cent for M-M gas , (3) 6+10 per 
cent for M-T, and (4) 25 ± 32 per cent for L\, \-M, respectively. 
At rJ o' X ' we found a similar difference in the normalization be- 
tween simulations and observations using either X-ray masses 
or weak lensing masses. 

We also derived the weak lensing mass based and X-ray 
mass based mass-observable relations at r Y x ' x without the LSS 
evolution correction (Figs.l9TfT0l>. For M-Yx, M-M gas , the con- 
clusions are relatively similar as for those with the LSS evolu- 
tion correction as the LSS evolution correction has the same 
impact on both variables. However, the LSS evolution cor- 
rection gives an opposite impact on Lboi and Fx (also M). 



Therefore the Lboi-Tx relation comparison between our sam- 
ple and the sample in Maughan et al. (2007) shows that the LSS 
evolution correction is required to account for the evolutionary 
effect of the scaling relations and the scatter is significantly re- 
duced using Yx measured at r Y x ' x (Fig. |C.5l ). 



6.2.3. All measurements within rj *' wl 

We followed a similar procedure outlined in §|4j to derive itera- 
tively the radius rsoo for each cluster, but using the mass derived 
from the weak lensing data. Hereafter we call it j*™" 1 . 

The procedure is initiated by deriving the mass profile 
M Yx (r) from the profile Y x (r) = M gas (r) * 7o. 2 -o.5r 5 oo using the 
MjL—Yx relation (see Table|2]i. The radial scale, within which 
the density contrast derived from the mass profile M Yx (r) is 
500, is called r test . We re-fit the M M -Y X relation but using 
the weak lensing masses and Yx parameters at r test for the 19 
LoCuSS clusters with the slope fixed to 0.568. This best fit 
is used to re-derive the mass profile M Yx (r) from the profile 
Yx(r), and M Yx (r) can then be used to re-determine the radial 
scale r test , within which the density contrast, derived from the 
mass profile M Yx (r), is 500. We iteratively re-fit the M wl -T x 
relation using the weak lensing masses and Y x parameters at 
r test , and re-determine the radial scale r test using M Yx (r) till the 
best fit of the M wl -ix relation does not vary any more for our 
sample. The radial scale, r™ w , is the final r test . The X-ray and 
weak lensing masses at r Yx ' wl are listed in Table. |C2l The X- 
ray mass and weak lensing mass based mass-observable re- 
lations at r Yx wl are shown in Figs. ICTHC21 and Table ED The 
normalization of the mass-observable relations using the X-ray 
masses and X-ray observables measured at r Y x ' wl is lower than 
the one from simulations by, i.e. (1) 20 + 3 per cent for M-Yx, 
(2) 16+4 per cent for M-M gas , (3) 18+5 per cent for M-T, and 
(4) 13 + 31 per cent for L\, i-M, respectively. The normaliza- 
tion of the mass-observable relations using the weak lensing 
masses and X-ray observables measured at r Y x ' wl is lower than 
the one from simulations by, i.e. (1) 17 + 8 per cent for M-Yx, 
(2) 14 + 8 per cent for M-M gas , (3) 14 + 10 per cent for M-T, 
and (4) 9 + 32 per cent for L\, \-M, respectively. 



6.2.4. All measurements within rj *' si 

We also derived rsoo, the radius by combing the T x (r) pro- 
file with the M-Yx relation from simulations in Nagai et al. 
(2007b). The M 5 oo-T x relation in Nagai et al. (2007b) can 
be used to derive the mass profile M Yx (r) from the profile 
Yx(r) = M gas (r) * 7o.2-o.5i 50 o- The radius, within which the 
density contrast derived from the mass profile M Yx (r) is 500, 
is the radial scale. Hereafter we call it r Yx si . The X -ray and 
weak lensing masses at r Y x ' wl are listed in Table. IC.3I The X- 
ray mass and weak lensing mass based mass-observable re- 
lations at rJ o' sl are shown in Figs. IC.3HC.4I and Table [5] The 
normalization of the mass-observable relations using the X-ray 
masses and X-ray observables measured at r™' sl is lower than 
the one from simulations by, i.e. (1) 15 + 3 per cent for M-Yx, 
(2) 16 + 5 per cent for M-M gas , (3) 8 + 5 per cent for M-T, and 
(4) 26 + 31 per cent for Lb<,i-M, respectively. The normaliza- 
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tion of the mass-observable relations using the weak lensing 
masses and X-ray observables measured at r^ sl is lower than 
the one from simulations by, i.e. (1) 15 + 8 per cent for M-Yx, 
(2) 21+8 per cent for M-M gas , (3) 6 ± 10 per cent for M-T, 
and (4) 25 ± 32 per cent for L^-M, respectively. 

We thus draw similar conclusions concerning the scatter 
and normalization in the mass-observable relations using the 
weak lensing masses at rj*' wl and r™ sl , respectively, as for 
those measured at rj^ x . As precisely shown above, the nor- 
malization of the M-Yx, M-Mg as , and M—T relations using X- 
ray mass estimates is lower than the values from simulations 
by ~ 20 per cent at 3<x significance. In good agreement, these 
relations based on weak lensing masses also show lower nor- 
malization by ~ 20 per cent than the values from simulations 
but at ~ 2cr significance. 

In Fig. IC.5I we also summarize the luminosity-Fx rela- 

Y y y w] Y si 

tions at , r,*' and r 5Q *'' . The upper panels show that the 
LSS evolution correction accounts for the segregation of the 
Luminosity versus Yx relations between cluster samples at dif- 
ferent redshifts. 

7. Discussion 

7.1. Scatter and bias in weak lensing mass estimates 

Since we observed a large scatter in weak lensing masses for 
the sample, we summarize the possible sources which can en- 
large scatter in weak lensing mass estimates as follows. 

(1) Accuracy of the faint galaxy shape measurements. 
Residuals in the removal of the point spread function (PSF) 
can masquerade as spurious shear signal, and/or introduce 
strong systematic errors, often referred to as B-models. Future 
LoCuSS weak lensing studies will use multiple shape measure- 
ment techniques to validate the fidelity of the PSF removal. 

(2) Size of FOV. Ideally the FOV should extend beyond the 
virialized region of each cluster, which for a 1O 15 M cluster at 
z^0.2 corresponds to r2oo-2.0Mpc, or >10'. Small FOV weak 
lensing observations typically incur systematic errors in clus- 
ter mass determination due to the mass-sheet degeneracy and 
the need to extrapolate the observed shear signal out to larger 
radii. Dahle's (2006) data are a good example of this; the NOT 
and UH 88in observations reach a maximum radius of just ~T. 
In contrast, the LoCuSS pilot study used wide field data from 
CFHT (Bardeau et al. 2007), and LoCuSS is using Suprime- 
CAM on Subaru to make weak lensing measurements (Okabe 
et al., in prep.) - both reach out to r>!5'. This will therefore 
not be a major concern for future LoCuSS papers. 

(3) Accuracy of the background galaxy redshifts. The weak 
lensing mass estimate depends on the assumed redshift distri- 
bution of background galaxies, including the extent to which 
the background galaxy catalog is contaminated by faint clus- 
ter members and foreground galaxies. Several color-selection 
techniques have been developed recently (Broadhurst et al. 
2005; Limousin et al. 2007) to deal with this issue. We will 
apply these methods within LoCuSS weak lensing analysis and 
explore refined techniques to measure reliable photmetric red- 
shifts to minimize contamination of the shear catalogs. 



(4) Projection of LSS. The scatter in the weak lensing 
masses could be slightly boosted by effects of projection due 
to LSS. For example, Hoekstra (2003) investigated effects of 
uncorrelated LSS on weak lensing mass measurements, show- 
ing that the accuracy of the measured mass of an NFW model 
does not improve once the outer radius for the fitting reaches 
the angular scale ~ 15'. Fortunately, the LoCuSS clusters are 
at redshift around 0.2, which places their virial radii at r<13'. 
A similar investigation of the impact of correlated LSS is also 
required. 

(5) Mass models. The weak lensing mass measurements 
using a tangential shear fitting are dependent on the assumed 
mass model, e.g. NFW model or SIS model. Usually an NFW 
model is considered to be more representative for the mass 
distribution in clusters than an SIS model. However, for some 
clusters, especially merging clusters, a NFW model is not fa- 
vored statistically, and the SIS model has only been ruled out at 
large significance in a small number of cases (e.g. Kneib et al., 
2003). LoCuSS weak lens modeling will test for the presence of 
substructures and attempt to fit multi-component models. The 
superb image quality of Subaru, and consequently higher num- 
ber density of background galaxies will also increase the prob- 
ability that we can discriminate between NFW and SIS models. 

(6) Effects of triaxiality and interior substructures. Clowe 
et al. (2004) investigated effects of cluster triaxiality and sub- 
structures on the mass measurements and found that the triaxi- 
ality is more dominant. The value of r2oo obtained by fitting the 
shear profile with an NFW model can vary by 10-15 per cent 
due to different projection, which corresponds to a mass (M200) 
change of ~ 40 per cent. King et al. (2001) show that the mass 
measurement can be affected by interior substructures by up to 
10 per cent. 

7.2. Cluster mass comparison between observations 
and simulations 

The cluster mass estimates using different approaches can be 
different as e.g. boost effects caused by merging systems might 
be different on the X-ray and weak lensing masses. For the 
purpose of precise cluster cosmology, one requires accurate 
knowledge of the systematics in X-ray cluster mass estimates 
and weak lensing cluster mass estimates and the difference be- 
tween cluster mass estimates from simulations and observa- 
tions. 

Despite the larger individual measurement uncertainties in 
weak lensing masses, the normalization agrees between mass- 
observable relations using X-ray mass estimates and weak lens- 
ing mass estimates. The mean of the weak lensing to X-ray 
mass ratios is about 1 with up to 51 per cent scatter (Fig.|7]and 
Fig.©. 

Therefore the X-ray and weak lensing approaches provide 
a very valuable tool to cross calibrate the normalization and 
to constrain the systematics of the mass-observable relations. 
The systematics in the calibration of cluster mass estimates can 
be reduced by using a sufficiently large sample of galaxy clus- 
ters with good statistic. The average of the weak lensing ver- 
sus X-ray mass discrepancy is about two times the errors of 
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the weak lensing masses (Fig. |7j. This indicates either that the 
true errors in the weak lensing masses might be two times as 
large as the value quoted or that the difference is a real dis- 
crepancy. Therefore, one really needs a sufficiently large sam- 
ple of galaxy clusters with both X-ray and weak lensing mass 
measured by large FOV observations to reduce the scatter in 
weak lensing mass estimates. This point will be more obvious 
with the progress of the LoCuSS project - the completion of 
the X-ray (XMM-Newton) and weak lensing (HST and Subaru) 
observations. 

At the radial scale of rj *j x , both the X-ray determined 
masses and the weak lensing determined masses systematically 
lie below the prediction from the mass-observable relations 
(i.e. M-Yx, M-Mg as and M-T) from simulations (Borgani et 
al. 2004; Kravtsov et al. 2006; Nagai et al. 2007b) by up to 
~ 20 per cent with ~ 3<x significance for the X-ray masses 
and ~ 2<t significance for the weak lensing masses (Figs. [5]- 
Figs. IC.ll4C.41 i. If one believes that the two independent ap- 
proaches, weak lensing and X-ray, give reliable mass estimates, 
more real physics has to be included in simulations to repro- 
duce the real galaxy clusters. This can be further verified by 
comparing the mass-observable relations between simulations 
and observations using a sufficiently large sample of galaxy 
clusters with both X-ray and weak lensing mass measurements. 
Since the difference in the normalization is only up to 20 per 
cent between mass-observable relations from observations and 
simulations, it is important to pursue a few per cent level sys- 
tematics in the mass estimates. Therefore one has to improve 
the statistic for both simulation and observational samples by 
using a large sample (e.g. ~ 100 galaxy clusters). And one also 
has to reduce the errors in observations by using both (1) deep 
X-ray observations (e.g. an accuracy of the ICM temperature 
measurements on the per cent level precision for 8-10 radial 
bins) performed with a large FOV (e.g. covering the area of 
more than rsoo), and (2) weak lensing observations performed 
with an instrument with a good PSF and a large FOV (e.g. 
covering the area of more than r2oo) together with a photo-z 
database to characterize the background source galaxy redshift 
distribution. 

7.3. Mass calibration and cr 8 

Based on (1) the number density of clusters with T > 6 keV 
(Henry 2004, Ikebe et al. 2002) and (2) three bins of the clus- 
ter density with certain velocity dispersion from Rines et al 
(2007), Evrard et al. (2007) found a high amplitude of cosmic 
power spectrum, S% = o"8(£2 m /0.3) ' 35 = 0.8 using the mass- 
temperature relation (M'-T') from simulations combined with 
the gas mass-temperature relation (M'—T') from ROS AT and 
ASCA observations in Mohr et al. (1999) also regarding mas- 
sive clusters in a similar temperature range. Our comparison 
shows that the normalization of the mass-temperature relation 
in Evrard et al. (1996) is higher than ours (M* QQ -T) by 40 per 
cent, and the normalization of the gas mass-temperature rela- 
tion in Mohr et al. (1999) is higher than ours (M gaS £oo~T) by 
8 1 per cent. The gas mass to mass ratio in Evrard et al. (2007) 
is M'/M' = (1.81/1.4) x Mg as /M. According to the gas mass 



to mass ratio, M' g JM' = 0.083 fc^Sg 5/3 , in Evrard et al. 

(2007), we obtained M gas /M = 0.083 x (1.4/1.81) h~^ 5 S & 5/3 . 
For our sample, we then found S% = 0.8571745^, where S'^ is 
0.8 for the sample in Evrard et al. (2007). We thus obtained 
S 8 = 0.69 for our sample. The X-ray mass calibration us- 
ing for our sample supports a lower amplitude of cosmic 
power spectrum on the basis of the comparison between obser- 
vations and simulations. 

To be conservative, we repeated the above calculation us- 
ing the mass-temperature and gas mass-temperature relations 
at the radius r J g sl , which is determined by our Yx profile and 
the mass-Tx relation from simulations in Nagai et al. (2007b). 
In this case, the normalization of the mass-temperature re- 
lation in Evrard et al. (1996) is higher than ours by 28 per 
cent, and the normalization of the gas mass-temperature re- 
lation in Mohr et al. (1999) is higher than ours by 58 per cent. 
Therefore the gas mass to mass ratio in Evrard et al. (2007) 
is M' g . d JM' = (1.58/1.28) M ga JM. According to the gas mass 

to mass ratio, M' gRS /M' = 0.083 h^ 5 S'^ 5/3 , in Evrard et al. 
(2007), we obtained M gas /M = 0.083x(1.28/1.58) h^ 5 S & 5 ' 3 . 
For our sample, we then found S& = 0.885 g = 0.70. 

Assuming weak lensing masses agree better with cluster 
masses from simulations, we also derived the S 8 parameter 
using the weak lensing mass based mass-temperature relation 
(M™Q Q -To.2-o.5r 500 ) with the slope fixed to 1.5, which normaliza- 
tion is 4.49 x 10 13 M o . It gives S 8 = 0.77, similar to the value 
obtained from simulations in Evrard et al. (2007). 

The above three experiments show the importance of the 
calibration of the scaling relations for the cluster cosmology. 
Depending on the mass calibration (e.g. via X-ray or weak 
lensing approach) the derived value of the amplitude of cos- 
mic power spectrum varies by 10 per cent between 0.69 and 
0.77 for our sample. 

8. Summary and conclusions 

We performed a systematic analysis to measure the X-ray quan- 
tities based on XMM-Newton observations for a sample of 37 
LoCuSS clusters. We investigated (1) the empirical similarity 
of the X-ray scaling relations, (2) the scatter in X-ray and weak 
lensing mass estimates, (3) the difference in the normalization 
between mass-observable relations from observations (using 
X-ray masses and weak lensing masses, respectively) and sim- 
ulations, and (4) the impact of the mass calibration on the de- 
termination of the amplitude of cosmic power spectrum. We 
summarize the main conclusions as follows. 

(i) The self-consistent X-ray scaling relations for our sample 
support an empirical similarity for massive galaxy clusters. The 
X-ray scaling relations of our sample show no evident evolu- 
tion compared to the nearby and more distant samples after the 
redshift evolution correction due to LSS growth. Specially, the 
X-ray M-T relations agree on 1 per cent level between our 
sample and the samples in Arnaud et al. (2005) and Vikhlinin 
et al. (2006a). 

The cluster cores (< 0.2r5oo) contribute up to 80 per cent to 
the bolometric X-ray luminosity. Using X-ray luminosity cor- 
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rected for the cool core and temperature excluding the cool 
core, leads both to less scatter (7-40 per cent) in the X-ray 
scaling relations and to better than 17 per cent agreement of 
the normalization between the sample and the non-CCC sub- 
sample. 

(ii) Using a Gaussian fit to the histogram of the weak lensing 
to X-ray M 500 ratios, the mean weak lensing to X-ray M 500 
ratios is about 1, indicating a good agreement between X-ray 
and weak lensing mass estimates for most clusters. Therefore 
the weak lensing approach can be a very valuable tool to cross 
calibrate the systematics in cluster mass estimates and the nor- 
malization of the mass-observable relations together with the 
X-ray approach. The comparison between weak lensing and 
X-ray mass estimates shows a large scatter of the individual 
results. The most robust way to characterize this scatter is by 
the width of the Gaussian fit to the distribution of deviations 
which is 31-51 per cent. The large scatter is mostly due to un- 
certainties in the weak lensing mass caused in part by the LSS 
into which the clusters are embedded. For the purpose of pre- 
cise cosmology, one requires a sufficiently large sample with 
both X-ray and weak lensing measurements to verify (1) the 
scatter in X-ray and weak lensing mass estimates on a few per 
cent level precision, and (2) the difference in the normalization 
between mass-observable relations between observations and 
simulations on per cent level accuracy. 

(iii) The scatter in the mass-observable relations is larger us- 
ing weak lensing masses than using X-ray masses by a factor 
of 2. However, the scatter using weak lensing masses neither 
depends on being classified as CCCs and non-CCCs nor de- 
pends on being classified as "single"s and "merger's. Using our 
method, we can define a radial scale r 5 ^ , which can not only 
reduce the scatter in the mass-observable scaling relations, but 
also provide consistency between weak lensing mass based and 
X-ray mass based mass-observable scaling relations. 

(iv) The X-ray mass estimates are about 20 per cent lower than 
the prediction of the mass-observable relations (i.e. M-Yx, M— 
M gas and M-T) from simulations with ~ 3cr significance. The 
comparison between the observed mass-observable relations 
(i.e. M-Y x , M-M gas and M-T at rj x ' x , rj x ' wl or rj x ' si ) and the 
relations from simulations (Borgani et al. 2004; Kravtsov et al. 
2006; Nagai et al. 2007b) shows that weak lensing mass esti- 
mates, similar to X-ray mass estimates, are about 20 per cent 
lower than the cluster masses predicted by the mass-observable 
relations from simulations, but with ~ 2<x significance. Since 
the X-ray and weak lensing approaches are independent, this 
might indicate that more complex physics is required in simu- 
lations to reproduce the real physics in galaxy cluster. 

(v) With a fixed slope to 1, the normalization of the best- 
fit power law for the M wl (rJ x ' x ) versus M x (rJ x ' x ) relation is 
1.09 + 0.08. On average, the weak lensing masses, M wl (rJ Q x ' x ), 
are higher than the X-ray masses, M x (rJ x ' x ), by 9 + 8 per cent 
according to the normalization. Following the arguments pre- 
sented in Churazov et al. (2007), the current comparison limits 



the non-thermal pressure support to 9 ± 8 per cent according to 
the ratio between the weak lensing and X-ray mass estimates. 

(vi) Depending on the mass calibration (e.g. via X-ray or weak 
lensing approaches), the derived value of the amplitude of cos- 
mic power spectrum varies by 10 per cent between 0.69 and 
0.77 for our sample. 

Since the weak lensing mass determination is expected to 
be unbiased on average for a large sample of clusters, it is also 
very important to compare the agreement between lensing and 
X-ray mass for the sample average. We find that the results are 
consistent with an agreement of both methods. But different ap- 
proaches to characterize the distribution of the deviation gives 
slightly different results. Taking these differences as a very con- 
servative measure of systematic uncertainties we conclude that 
the two methods of mass determination agree within an uncer- 
tainty of less than 20 per cent for the sample average. 

It is clear from this result that the latter conclusion, which is 
very important for the cosmological application of clusters, will 
benefit very much from an increase of the sample size, in addi- 
tion to the currently happening improvements of lensing mass 
determinations. The present paper therefore highlights the cru- 
cial importance of the completion of the entire LoCuSS sample, 
which is planned to comprise in the end about 100 clusters. 
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Table 1. Power law, Y = Yq X 7 , parameterized X-ray scaling relations. The mean of the scatter is determined by a Gaussian 
model fit to the histogram of the scatter. 
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Here we provide the mean scatter derived by a Gaussian fit to the histogram, in which double values are from the histogram with a bi-modal 
Gaussian distribution. In Zhang et al. (2006, 2007), we only presented the average scatter due to the limited data points. 
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Table 2. Power law, Y = Yq X y , parameterized mass-observable scaling relations using weak lensing masses (MfL) for the 19 
LoCuSS clusters having weak lensing masses. 
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In Col. Sample, "merger" (which denotes "primary with small secondary", "elliptical", "off-center" and "complex") and "single" are clas- 
sified as described in Jones & Forman (1992). Here we provide the mean scatter derived by a Gaussian fit to the histogram, in which double 
values are from the histogram with a bi-modal Gaussian distribution. 

Table 3. See the caption in Table [2] except that the weak lensing masses, X-ray masses and Yx parameters are determined within 
the radii (rJ x ' X ) determined by the method described in § [4] 
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Table 4. See the caption in Table [2] except that the weak lensing masses, X-ray masses and Fx parameters are determined within 
the radii ('V ) x ' wl ) determined by the method described in § 16.2.31 
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Table 5. See the caption in Table [2] except that the weak lensing masses, X-ray masses and Fx parameters are determined within 
the radii (r 5 £ sl ) determined by the method described in § 16.2.41 
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M ,as.500 E ( z ) M gas.500 E ( z ) ( M S ) 

Fig. 1. Mass-ix (upper panels) and mass-gas mass (lower panels) relations using X-ray masses (left panels) and weak lensing 
masses (right panels). The CCCs are marked by triangles in the left panels. The best-fit power law in the right panels (solid, 
using weak lensing masses; dashed, using X-ray masses) with their slopes fixed to the values as found in simulations in order 
to compare the normalization (Table|2]i. The clusters appearing "primary with small secondary" or "off-center" morphology are 
in green, appearing "elliptical" or "complex" morphology are in blue, and appearing "single" morphology are in red using the 
classification in Jones & Forman (1992). 
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Fig. 2. Mass-temperature (upper panels) and luminosity-mass (lower panels) relations using X-ray masses (left panels) and weak 
lensing masses (right panels). The CCCs are marked by triangles in the left panels. The best-fit power law in the right panels 
(solid, using weak lensing masses; dashed, using X-ray masses) with their slopes fixed to the values as found in simulations in 
order to compare the normalization (Table|2]i. The colors have the same meaning as those in Fig.Q] 




Fig. 3. Left: Gas mass-temperature relation. Right: Luminosity- Yx relation. The CCCs are marked by triangles. The colors have 
the same meaning as those in Fig.Q] 




Fig. 4. X-ray luminosity in the 0.1-2.4 keV band versus temperature (left) and bolometric X-ray luminosity versus temperature 
(right). The CCCs are marked by triangles. The colors have the same meaning as those in Fig. Q] 




V^oo™)^) (Me) M gas (r 500 Y «' x )E(z) (M ) 



Fig. 5. See caption in Fig. Q] except that the weak lensing masses, X-ray masses and Y x parameters are determined within the 
radii (r™ ) determined by the method described in §|4](also see TableO. 
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Fig. 6. See caption in Fig. [2] except that the weak lensing masses, X-ray masses and Yx parameters are determined within the 
radii (r™ ) determined by the method described in §|4](also see Table[3j. 
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Fig. 7. Upper left: M" ' versus M* (a best-fit slope of 1.61 + 0.39). Upper right: M^ Q versus M^* (a best-fit slope of 
1.57 ± 0.38). Mj * is derived from Yx with the slope fixed to 0.568 (the prediction from simulations in Nagai et al. 2007b) fitting 
M-Yx relation for our sample. Lower panels: Normalized differential cluster number count as a function of the M™ ' to M* 00 
ratio and M^J () to ratio, respectively, with their Gaussian fits (gray, solid curves denote the fits with x-axis in linear space, 
and dashed curves in logarithmic space). The lower parts of the lower panels show the M™ ' () versus Mf 00 discrepancy and M^J 
versus discrepancy, respectively, normalized to the error bar of M?L. The colors have the same meaning as those in Fig.Q] 
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Fig. 9. See caption in Fig. [5]except that no LSS evolution correction is applied in the mass-observable relations. 
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Fig. 10. See caption in Fig. [6]except that no LSS evolution correction is applied in the mass-observable relations. 
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Appendix A: X-ray mass modeling 

The X-ray mass modeling described here helps to understand 
the measurement uncertainties of the X-ray quantities used for 
studies of the scaling relations. 

We show the ROSAT luminosity versus redshift distribu- 
tion for the LoCuSS sample in Fig. lA.ll and the observational 
information of the XMM-Newton data for the 37 LocuSS cluters 
in Table lA.fl The primary parameters of all 37 LoCuSS galaxy 
clusters are given in Table I A. 2l 



A.1. Density contrast 

To determine the global cluster parameters, one needs a fiducial 
outer radius that was defined as follows. The mean cluster den- 
sity contrast is the average density with respect to the critical 
density, 



A«r) 



3M(< r) 
4-nr 3 p c (z) 



(A.l) 



The critical density at redshift z is p c (z) = p c oE 2 (z), where 
E 2 (z) = O m (l + z) 3 + Q A + (1 - Q m - Qa)(1 + zf. r A is the 
radius within which the density contrast is A. Ma is the total 
mass within r&. For A = 500, rsoo is the radius within which 
the density contrast is 500 and M500 is the total mass within 
rsoo- 

A.2. Temperature profiles 

XMM-Newton (as well as Chandra), in contrast to earlier tele- 
scopes, has a less energy-dependent, smaller PSF, more reliable 
to study cluster temperature profiles. We de-projected the spec- 
tra (see the method in Zhang et al. 2007) and performed the 
spectral fitting in XSPEC to obtain the radial profiles of the 
temperature and metallicity. In Figs. |A.2jlA.4| in the electronic 
edition of the Journal (together with Fig. 2 in Zhang et al. 2006, 
and Fig. 1 and Figs. B. 1-B.5 in Zhang et al. 2007 for the whole 
sample) we show the radial temperature profiles. The tempera- 
ture profiles are approximated by 

T(r) = P 3 exp[-(r - Pif/P 2 ] + P 6 (l + r 2 /P 2 4 y P5 + P 7 , (A.2) 

motivated by the qualitative appearance of the observed tem- 
perature distributions of clusters of galaxies. P,-, i = 1, 7, are 
used for parameterization. 

A CCC often shows both a peaked surface brightness pro- 
file and a steep temperature drop towards the cluster center. Our 
sample does not show pronounced CCC and non-CCC bimodal 
distribution. For studies of the segregation of the X-ray scaling 
relations we empirically define the CCC as the cluster whose 
temperature drop is greater than 30 per cent of the peak tem- 
perature towards the cluster center. Nine clusters are therefore 
CCCs in this sample. There might be more CCCs in the sample 
if the sample would be constructed using the higher resolution 
temperature measurements. 

The X-ray surface brightness isophotes are also sensitive 
to cluster populations. We therefore also used the classifica- 
tion of the cluster morphologies in Jones & Forman (1992) and 
divided the sample to 5 types (Table |B~TT i: "single", "primary 



with a small secondary", "elliptical", "off-center" and "com- 
plex" based on the smoothness of its X-ray surface brightness 
isophotes in the 0.7-2 keV band. Therefore "single" approx- 
imately denotes relaxed clusters, "primary with a small sec- 
ondary", "off-center", "elliptical" and "complex" denote merg- 
ing clusters. 

About 80 per cent CCCs are relaxed clusters ("single" de- 
fined in §[6] using surface brightness isophotes). As we show 
later, these 9 CCCs often show shorter central cooling time, 
larger cooling radii (scaled to rsoo) and lower central entropies 
compared to the non-CCCs in the sample (Table |B~TT >. 

A.3. Surface brightness profiles 

A /? model (e.g. Cavaliere & Fusco-Femiano 1976; Jones & 
Forman 1984) is often used to describe electron density pro- 
files in clusters. To obtain an acceptable fit for all clusters in the 
sample, we adopt a double-/? model of the electron number den- 
sity profile, « e (r) = n e0 i(l + r 2 1 r 2 ^' 2 + « e02 (l + r 2 /r 2 ^ 2 , 
where n e o = w e oi + n e02 is the central electron number density 
(see Table IB. il l. /? the slope parameter, and r c \ and r c2 the core 
radii of the inner and outer components, respectively. The /? 
value is derived by fitting the outer component, and it is used 
for the inner component since the inner component is insensi- 
tive to the slope. The soft band (e.g. 0.5-2 keV for 9 REFLEX- 
DXL clusters and 0.7-2 keV for the remaining 28 clusters) X- 
ray surface brightness profile model S(R), in which R is the 
projected radius, is linked to the radial profile of the ICM elec- 
tron number density n e (r) and emissivity function as an integral 
performed along the line-of-sight for galaxy clusters. For mas- 
sive clusters (i.e. T > 3.5 keV) the cooling function can be 
treated as a temperature independent quantity in the soft band, 
and the surface brightness profile is thus described as follows, 



S X (R) <x 



«p« ( 

R 



dl. 



(A.3) 



For the whole sample, the best-fit "mekal" model of the spectra 
limited to 0.2r5oo < R < 0.5rsoo is used to convert the surface 
brightness profile to derive the density distribution. We fitted 
the observed surface brightness profile by this integral con- 
volved with the XMM-Newton PSF matrices (Figs. |A.5llA.7| in 
the electronic edition of the Journal, Fig. 3 in Zhang et al. 2006, 
and Fig. 1 and Figs. B.1-B.5 in Zhang et al. 2007 for the whole 
sample) and obtained the parameters of the double-/? model of 
the electron density profile. The fit was performed within the 
truncation radius (r t , see Table lA.2"b corresponding to a S /N of 
3 of the observational surface brightness profile. The truncation 
radii r t are often larger than rsoo- The cluster cores, referring to 
the core of the inner component of the double-/? model, span a 
broad range up to 0.2 rsoo- 

A. 4. Mass distribution 

The ICM can be used to trace the cluster potential. We assume 
that, (1) the ICM is in hydrostatic equilibrium within the grav- 
itational potential dominated by DM, (2) the DM distribution 
is spherically symmetric, and (3) the kinematic term is approx- 
imated to zero in hydrostatic equilibrium. The cluster mass is 
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then calculated from the X-ray measured ICM density and tem- 
perature distributions by, 

1 d[n e (r)T(r)] _ GM(< r) 

7\ 1 _ 2 ' (A- 4 .) 

fim p n e {r) dr r L 

where fi = 0.62 is the mean molecular weight per hydrogen 
atom. The Boltzmann constant is included in T as the temper- 
ature is in units of keV. Following the method in Zhang et al. 
(2007), we use a set of input parameters of the approximation 
functions, in which n e oi, r c i (i = 1,2) represent the elec- 
tron number density profile n e (r) and P\ (i — 1, 7) represent 
the temperature profile T(r), respectively, to compute the mean 
cluster mass. The mass uncertainties are propagated using the 
uncertainties of the electron number density and temperature 
measurements by Monte Carlo simulations as follows (see also 
Zhang et al. 2007). For each cluster, we simulated the electron 
density and temperature profiles of a sample of 100 clusters 
using observed electron density and temperature profiles and 
their errors. The mass profiles and other properties of the 100 
simulated clusters were calculated to derive the errors. 

The derived mass profile was used to calculate the den- 
sity contrast profile to derive M500 and rsoo (see Table IB. fl . In 
Figs. lA.8l4A.10l in the electronic edition of the Journal (together 
with Fig. 5 in Zhang et al. 2006, and Fig. 1 and Figs. B.1-B.5 
in Zhang et al. 2007 for the whole sample), we present the ob- 
served mass profiles. 



A. 5. Gas mass fraction distribution 

The gas mass fraction is an important parameter for cluster 
physics, e.g. heating and cooling processes, and cosmological 
applications using galaxy clusters (e.g. Vikhlinin et al. 2002; 
Allen et al. 2004). The gas mass fraction distribution is defined 
as /gas(< r) = M gas (< r)/M{< r). The gas mass, total mass 
and gas mass fraction at are given in Table IB. ll We obtain 
an average gas mass fraction of 0.121 + 0.004 at for the 
sample. 

As observed, the gas mass fraction at r^oo is generally lower 
than the universal baryon fraction, ft, = Q.\,/Q. m = 0.176 + 
0.019, based on the WMAP three year results in Spergel et 
al. (2006). This is because the baryons in galaxy clusters re- 
side mostly in hot gas together with a small fraction of stars as 
implied in observations (0.003-0.022 in groups, Gastaldello et 
al. 2007) and simulations (15 per cent of / gas , e.g. Kravtsov et 
al. 2005). In principle, Q m can be determined from the baryon 
fraction, /j, = / gas + / ga i, in which a contribution from stars in 
galaxies is given by / ga i = 0.02 + O.Olh^ (White et al. 1993). 
The average gas mass fraction at /-500 for the sample supports a 
low matter density Universe as also shown in e.g. Allen et al. 
(2002), Ettori et al. (2003) and Vikhlinin et al. (2003). 
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+32 


: 07 : 


47.4 


0093031001 


2003-08-29 


Thin 


Thin 


EFF 


2.7 


2.9 


0.9 


Abell2390 


21 : 


: 53 


: 37.1 


+ 17 


: 41 : 


46.4 


0111270101 


2001-06-19 


Thin 


Thin 


FF 


10.3 


10.0 


8.8 


Abell2667 


23 : 


: 51 


: 39.2 


-26 


: 05 : 


03.5 


0148990101 


2003-06-21 


Medium 


Medium 


FF 


22.2 


22.9 


14.1 


Z7160 


14 : 


: 57 


: 15.2 


+22 


: 20 : 


31.2 


0108670201 


2002-08-03 


Medium 


Medium 


FF 


31.2 


31.6 


26.3 



The cluster center is in sky coordinates in epoch J2000. The MOS data are in EE mode except that the MOS1 data of Abelll835 are in 
window mode, which cannot be used for this work. 
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Table A.2. Primary parameters. Column (1): cluster name; Col. (2): optical redshift (e.g. Bohringer et al. 2004; Smith et al. 
2005; Bardeau et al. 2007); Col. (3): hydrogen column density (Dickey & Lockman 1990); Col. (4): truncation radius corre- 
sponding to a S IN of 3 of the observational surface brightness profile; Col. (5): volume averaged radial temperature profile 
of 0.2-0.5r 5 oo; Cols. (6-7): Spectroscopic temperature and metallicity in the annulus of 0.2-0.5r 5( x). Cols. (8-10): bolometric 
luminosity including the < 0.2r 5 oo region, excluding the < 0.2rg o region, and corrected for the < 0.2r 5 oo region, respectively. 



Name 




N a 


n 


T<a-0,5r m 


rj-,spec 
0.2-0.5,™ 


Z 0.2-0.5,™ 


j incc 
Ho] 


j exec 
Hoi 


j corr 
Hoi 






10 22 cm" 2 


t 


keV 


keV 


z Q 


10 45 erg/s 


10 45 erg/s 


10 45 erg/s 


RXCJ0043. 4-2037 


0.292 


0.0154 


7.0 


6.8 ± 0.4 


7.0 


±0.4 


0.17 ±0.06 


1.05 ±0.10 


0.59 ± 0.07 


0.84 ±0.10 


RXCJ0232.2-4420 


0.284 


0.0250 


6.8 


6.9 ± 0.3 


6.6 


±0.3 


0.23 ± 0.07 


1.85 ±0.14 


0.75 ± 0.08 


1.10 + 0.14 


RXCJ0307.0-2840 


0.258 


0.0136 


4.7 


6.9 ± 0.4 


7.1 


±0.4 


0.26 ± 0.08 


1.31 ±0.12 


0.61 ±0.08 


0.88 + 0.12 


RXCJ05 16.7-5430 


0.294 


0.0686 


6.4 


6.7 ± 0.5 


6.7 


±0.5 


0.25 ±0.10 


1.09 ±0.15 


0.80 ±0.11 


1.00 + 0.15 


RXCJ0528.9-3927 


0.284 


0.0212 


5.8 


6.6 ± 0.5 


6.6 


±0.5 


0.20 ± 0.09 


1.45 ±0.15 


0.88 ±0.10 


1.10 + 0.15 


RXCJ0532.9-3701 


0.275 


0.0290 


7.0 


8.2 ± 0.6 


7.7 


±0.6 


0.28 ± 0.09 


1.27 ±0.12 


0.55 ± 0.07 


0.81 ±0.12 


RXCJ0547.6-3152 


0.148 


0.0205 


7.7 


6.1+0.2 


6.0 


±0.2 


0.22 ± 0.03 


0.72 ± 0.06 


0.39 ± 0.04 


0.55 ± 0.06 


RXCJ0645.4-5413 


0.164 


0.0651 


8.1 


8.1+0.3 


7.6 


±0.3 


0.11 ±0.04 


1.78 ±0.16 


1.01+0.12 


1.30 + 0.16 


RXCJ0658.5-5556 


0.296 


0.0653 


8.9 


11.1+0.4 


10.7 + 0.4 


0.14 ±0.04 


4.95 ± 0.24 


2.99 + 0.18 


4.21 ±0.24 


RXCJ0945.4-0839 


0.153 


0.0359 


7.6 


5.7 ± 0.5 


5.3 


±0.5 


0.11 ±0.10 


0.42 ± 0.06 


0.24 ± 0.04 


0.33 ± 0.06 


RXCJ0958.3-1103 


0.167 


0.0540 


7.4 


6.1+0.3 


5.8 


±0.3 


0.27 ± 0.06 


0.91+0.09 


0.42 ± 0.06 


0.57 ± 0.09 


RXCJ2 129.6+0005 


0.235 


0.0428 


6.8 


6.4 ± 0.2 


6.3 


±0.2 


0.27 ± 0.03 


1.43 ±0.09 


0.58 ± 0.06 


0.79 ± 0.09 


RXCJ2218.6-3853 


0.141 


0.0138 


7.3 


6.7 ± 0.2 


6.2 


±0.2 


0.27 ± 0.05 


0.66 ± 0.05 


0.34 ± 0.04 


0.48 ± 0.05 


RXCJ2234.5-3744 


0.151 


0.0122 


7.9 


7.8 + 0.4 


7.5 


±0.4 


0.22 ± 0.08 


1.09 ±0.11 


0.57 ± 0.07 


0.80 + 0.11 


RXCJ2308.3-0211 


0.297 


0.0445 


7.3 


7.7 ± 0.7 


7.6 


±0.7 


0.41 ±0.16 


1.21 ±0.13 


0.56 ± 0.08 


0.80 ±0.13 


RXCJ2337.6+0016 


0.275 


0.0382 


7.0 


8.0 ± 0.5 


7.5 


±0.5 


0.19 ±0.07 


1.00 ±0.09 


0.54 ± 0.06 


0.77 ± 0.09 


Abell68 


0.255 


0.0493 


6.1 


7.7 ± 0.3 


7.3 


±0.3 


0.12 ±0.05 


1.14 ± 0.10 


0.61 ±0.07 


0.86 + 0.10 


Abellll5 


0.197 


0.0543 


8.0 


6.4 + 0.1 


6.2 


±0.1 


0.20 ± 0.03 


1.43 ±0.11 


1.05 ±0.09 


1.16 + 0.11 


Abell209 


0.209 


0.0164 


7.4 


7.1+0.3 


7.1 


±0.3 


0.22 ± 0.04 


1.33 + 0.11 


0.84 ± 0.08 


1.10 + 0.11 


Abell267 


0.230 


0.0280 


5.3 


6.5 ± 0.4 


6.2 


±0.4 


0.25 ± 0.08 


0.77 ± 0.07 


0.36 ± 0.05 


0.52 ± 0.07 


Abell383 


0.187 


0.0392 


7.7 


5.3 + 0.2 


4.7 


±0.2 


0.18 ±0.06 


0.81+0.05 


0.27 ± 0.04 


0.36 ± 0.05 


Abell773 


0.217 


0.0144 


7.7 


8.1+0.4 


8.3 


±0.4 


0.30 ± 0.06 


2.09 ±0.16 


1.06 + 0.11 


1.53+0.16 


Abell781 


0.298 


0.0194 


6.5 


6.4 ± 0.5 


6.5 


±0.5 


0.18 ±0.10 


0.63 ±0.10 


0.52 ± 0.09 


0.60 + 0.10 


Abell901 


0.163 


0.0499 


4.6 


3.6 ± 0.2 


3.2 


±0.2 


0.12 ± 0.09 


1.60 ± 0.05 


0.35 ± 0.03 


0.57 ± 0.05 


Abell963 


0.206 


0.0140 


6.1 


6.3 ± 0.2 


6.5 


±0.2 


0.28 ± 0.05 


1.14 ±0.09 


0.55 ± 0.06 


0.74 ± 0.09 


Abelll413 


0.143 


0.0219 


8.4 


6.8 + 0.1 


6.6 


±0.1 


0.27 ± 0.03 


1.38 ±0.08 


0.68 ± 0.06 


0.93 ± 0.08 


Abelll689 


0.184 


0.0182 


8.0 


8.4 ± 0.2 


8.5 


±0.2 


0.26 ± 0.04 


2.84 ±0.10 


0.93 ± 0.06 


1.44 + 0.10 


Abelll758 


0.280 


0.0106 


6.0 


7.6 ± 0.2 


7.9 


±0.2 


0.19 ±0.04 


1.05 ±0.07 


0.66 ± 0.05 


0.94 ± 0.07 


Abelll763 


0.228 


0.0936 


7.5 


6.3 ± 0.3 


5.8 


±0.3 


0.20 ± 0.05 


1.65 ±0.15 


1.14 + 0.11 


1.41 ±0.15 


Abelll835 


0.253 


0.0232 


7.0 


8.0 ± 0.3 


8.4 


±0.3 


0.23 ± 0.04 


5.32 ±0.17 


1.51 ±0.10 


2.25 ±0.17 


Abelll914 


0.171 


0.0095 


8.2 


8.6 + 0.3 


8.8 


±0.3 


0.19 ±0.05 


2.17 ±0.11 


0.62 ± 0.06 


0.97 ±0.11 


Abell2204 


0.152 


0.0567 


9.0 


8.0 ± 0.2 


7.6 


±0.2 


0.27 ± 0.03 


3.39 ±0.12 


1.09 + 0.08 


1.52 ±0.12 


Abell2218 


0.176 


0.0324 


7.7 


7.4 ± 0.3 


6.6 


±0.3 


0.21 ±0.04 


1.11 ±0.08 


0.68 ± 0.06 


0.90 ± 0.08 


Abell2261 


0.224 


0.0328 


8.1 


7.0 ± 0.6 


6.6 


±0.6 


0.33 ±0.18 


1.39 ±0.22 


0.67 + 0.13 


0.91 ±0.22 


Abell2390 


0.233 


0.0680 


7.7 


10.6 + 0.6 


11.6 + 0.6 


0.24 ± 0.06 


4.09 ± 0.27 


1.95 ±0.17 


2.66 ± 0.27 


Abell2667 


0.230 


0.0165 


7.6 


7.3 + 0.3 


7.0 


±0.3 


0.29 ± 0.04 


2.16 ± 0.13 


0.83 ± 0.08 


1.14 + 0.13 


Z7160 


0.258 


0.0316 


8.0 


5.0 + 0.1 


4.7 


±0.1 


0.32 ± 0.04 


1.70 ±0.08 


0.50 ± 0.05 


0.67 ± 0.08 



The global temperature and luminosity are uniformly calculated as described in Appendix lB.2l and Appendix lB.4l for all 37 LoCuSS clusters. 
In Zhang et al. (2006, 2007), the temperature and luminosity are defined differently. 
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0) 




Fig.A.l. ROSAT luminosity versus redshift for the LoCuSS 
sample. The crosses denote the 37 LoCuSS clusters observed 
by XMM-Newton. 

Appendix B: Scaled profiles of the X-ray 
properties: a structural similarity 

Simulations (e.g. Navarro et al. 1997, 2004) suggest a self- 
similar structure for galaxy clusters within hierarchical struc- 
ture formation scenarios. The scaled profiles of the X-ray prop- 
erties and their scatter can be used to quantify the structural 
variations, to probe the regularity of galaxy clusters, and thus 
to understand the scatter in the X-ray scaling relations, which 
is of prime importance for cosmological applications involving 
clusters of galaxies. 

Because the observational truncation radii (r t ) in the sur- 
face brightness profiles are above rsoo but below r2oo for most 
clusters in the sample, we use r^oo for radial scaling. 

The following redshift evolution corrections due to LSS 
growth are usually used to account for the dependence on the 
evolution of the cosmological parameters, 

S x -E- 3 (z)Kf(T), 

S ■ E A l\z) oc f(T), 

P ■ E-\z) oc f(T), 

L-E-\z)ocf(T), 

M ■ E(z) oc f(T), 

M gas • E(z) oc f(T\ 

Y x ■ E(z) oc f(T), 

r ■ E(z) oc f(T), 

where Yx = M gas • T is an X-ray analog for the integrated 
SZ flux (e.g. Kravtsov et al. 2006; also SZ see Sunyaev & 
Zeldovich 1972). The function f(T) denotes the best-fit power 
law parameterization. A c ,/A C; o factors dropped out compared 
to the redshift evolution corrections used in Zhang et al. (2006, 
2007) for convenience, and for easier comparison between our 
sample and published samples. 



3 A c , z = 18/r 2 + 82(n m . z - 1) - 39(n ra , z - l) 2 is the analytic ap- 
proximation derived from the top-hat spherical collapse model for a 
flat Universe (Bryan & Norman 1998) and Si mz the cosmic density 
parameter at redshift z. 



B.1. Scaled temperature profiles 

We scaled the radial temperature profiles by 7o.2-o.5r 50 o an d r 500 
(Fig. IB. 11 1. Within 0.2r5oo, we observed a temperature drop to 
at least 70 per cent of the maximum value towards the clus- 
ter center for 9 clusters (defined as CCCs). Previous obser- 
vations have shown that temperature measurements on scales 
below 0.2rgoo tend to show peculiarities linked to the cluster 
dynamical history (e.g. Finoguenov et al. 2005; Zhang et al. 
2005a, 2005b). For example, the temperatures of merging clus- 
ters can be boosted (Smith et al. 2005, Poole et al. 2006, 2007). 
However, the boosting affects mainly the cluster cores. Using 
global temperatures excluding the < O.2/-500 region can thus (1) 
tighten the scaled profiles of the X-ray properties, (2) reduce 
the scatter in the X-ray scaling relations, and (3) improve agree- 
ment between normalization of the X-ray scaling relations for 
the CCCs and non-CCCs, respectively. 

We used the 1-cr boundary of the scaled radial temper- 
ature profiles to compute the confidence zone for the mean 
temperature profile. The average temperature profile for the 
whole sample gives T(r) oc r I3±0 04 5 an( j f or the CCC sub- 



sample T(r) 



„0.38±0.04 



respectively, in the r < 0.2r5oo region. 



Because the temperature profiles are not PSF corrected, the sig- 
nificance of the cool cores could be more pronounced consid- 
ering the PSF effects. The temperature behavior in the cluster 
cores for the CCC subsample is very similar to the behavior 
for the nearby CCC sample in Sanderson et al. (2006) giving 
T(r) oc r 0A based on Chandra observations. In the outskirts 
(0.2r5oo < r < rsoo), the whole sample shows a structurally 



similar behavior giving T(r) > 



-0.28±0.19 



with 30 per cent scat- 



ter. The average temperature profile decreases down to 80 per 
cent of the maximum value with intrinsic scatter of ~ 30 per 
cent at ~ r^oo- 

Studies of the cluster temperature distributions indicate a 
similar universal temperature profile with a significant decline 
beyond cluster cores in observations (e.g. Markevitch et al. 
1998; De Grandi & Molendi 2002; Vikhlinin et al. 2005; Zhang 
et al. 2004, 2005a, 2006, 2007; Pratt et al. 2007) and simula- 
tions (e.g. Borgani et al. 2004). 

Accurate temperature distributions (15 per cent uncertain- 
ties) are measured up to 0.5r500 for our sample. We found 
steepening of the average temperature distribution using the pa- 
rameterization of the temperature distributions, giving slopes 
from 0.04 + 0.03 at 0.2r 50 o to 0.16 + 0.03 at 0.5r 5O0 (the er- 
rors for the slope are the statistical formal errors). The slope 
should have been calculated from the observational data points 
for the temperature profile as done for the density profile based 
on a temperature profile better resolved spatially. With deep 
Chandra observations, Vikhlinin et al. (2005) found more sig- 
nificant steepening of the temperature distributions in the out- 
skirts. 

B.2. Global temperature 

As indicated by the scaled temperature profiles for this sam- 
ple and also explained in Zhang et al. (2007), the temperature 
measurements limited to 0.2-0.5r5oo are less affected by the 
cool cores and are thus used as cluster global temperatures. We 



32 



Zhang et al.: LoCuSS: comparison of observed X-ray and lensing galaxy cluster scaling relations with simulations 




r (arcmin) 



r (arcmin) 



RXCJ0945. 4-0839 




r (arcmin) 



r (arcmin) 




r (arcmin) 



r (arcmin) 



Fig.A.2. De-projected temperature profiles approximated by the parameterization T(r) = P3 exp[-(r - P\) /P2] + P6(X + 
r 2 /P?)~ Ps + Pi crossing all the data points (solid). The vertical line denotes O.lrsoo. 



present the spectroscopic global temperature (To2^05r ) anc ^ ^ es °^ tne structural similarity of the X-ray properties and the 

metallicity (Z s ^l Q5r5 J in the annulus of 0.2-0.5r 5 oo together scaling relations. We present the spectral measurements, tem- 

wifh the volume average of the radial temperature profile lim- perature versus metallicity, from all radial bins in Fig. |B3J and 

ited to the radial range of 0.2-0.5r 500 (7o.2-o.5r 500 ) in Table[0 observed no obvious correlation between cluster temperature 

7o.2-o.5r 500 is used as the cluster global temperature for stud- and metallicity though the mass range is relatively small. 
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Fig. A.3. See caption in Fig. lA.2 



B.3. Scaled surface brightness profiles T 1A scaled surface brightness profiles in Fig. IB.2I We found a 

r ca , , less scattered structurally similar behavior at R > 0.2r5oo for 

The standard self-similar model gives Sx x L nidi oc M m , , , <■ . ■ £1 , . cl 

b JR e g JS the scaled surface brightness prohles compared to the profiles 

M~ 5 ' oc T as M oc M gas oc T . The scaling becomes scaled by 5x °c 5 
steeper (5x « r 1 ) using the empirical relations (M oc r 15 
in e.g. Arnaud et al. 2005, M gas oc r 1 8 in e.g. Mohr et al. 1999; 
Vikhlinin et al. 1999). We present the empirical scaling Sx x 
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Fig. A.4. See caption in Fig. IA.2l 
B.4. Global luminosities 

The bolometric X-ray luminosity (here we use the 0.01- 
100 keV band) is given by Lboi 00 f A.{r)n^.{r)dV , practically 
an integral of the X-ray surface brightness to 2.5r5oo (^™] C m 
Tabl dA.2l >. The value 2.5rsoo is used because the virial radii 
are about 2.2-2.6r 5 oo using our mass modeling. The luminos- 
ity varies by less than 10 per cent when setting the truncation 
radius to values between r 5 oo and 2.5r5oo- The best-fit "mekal" 
model for the spectra extracted in the annulus limited to 0.2- 
0.5r5oo is used to convert the count rate to real flux. Using dif- 
ferent converting "mekal" models, e.g. different temperature 
(7-0.2-0.5^0 or 7 , oT- C o.5r 500 or r o P 2-o.5r 500 or T <2/ir) and metallicity 

( Z ai-0.5r 500 or Z o P 2 e -o.5,- 50 „ OT Z <ik) causes U P to 20 P er cent de- 
ference in the X-ray luminosity, in which the difference is more 

significant for pronounced CCCs. 

The core radii populate a broad range of values up to 
0.2 rsoo. The X-ray luminosity is sensitive to the presence of 
a cool core. It can thus be used to probe the evolution of 
a cool cores. We present the bolometric luminosities includ- 
ing and excluding the R < 0.2rsoo region {l^f and L^f) in 
Table IA.2I The fractions of the X-ray luminosity attributed to 



the R < 0.2r5oo region span a broad range up to ~ 80 per 
cent. This introduces large uncertainties in the use of the lu- 
minosity LjJ 1 ^ as a mass indicator. The X-ray luminosity ex- 
cluding the < 0.2r 5 oo region (L^ c ) is much less dependent on 
the presence of the cool core. However, to better re-produce 
the cluster luminosity and the normalization of the luminosity- 
temperature/mass relation, we use the X-ray luminosity cor- 
rected for the cool core by assuming S x(R) — Sx(0.2r5oo) for 
the R < 0.2r m region (Lg™ in Table Using for the 
X-ray scaling relations we obtain reduced scatter and consistent 
normalization for the CCCs and non-CCCs as shown later. The 
luminosity is lower (by up to 10 per cent) assuming a constant 
luminosity in the cluster core instead of a (5 model, especially 
for the pronounced CCCs. However, the luminosity corrected 
for the cool core by the /3 model still introduces relatively sig- 
nificant scatter dominated by the CCCs due to the correlation 
between core radius and slope in the f3 model. 



4 This assumption was used only in the calculation for LJj°", not for 
any other quantity. 
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Fig. A.5. Surface brightness profiles parameterized by double-/? model. 



B.5. Gas profiles in the outskirts 

The generally adopted j3 model (ft = 2/3) gives n e oc r~ 2 . 
However, Vikhlinin et al. (1999) found a mild trend for /? 
to increase as a function of cluster temperature, which gives 
P ~ 0.80 and « e oc r~ 2A for clusters around 10 keV. Bahcall 
(1999) also found that the electron number density scales as 



n e oc r~ 2 4 at large radii. Zhang et al. (2006, 2007) confirmed 
their conclusion that n e oc r~ 2 A1 for the REFLEX-DXL clus- 
ters and « e oc r - 2 2±0 i f or the pilot LoCuSS sample. Due to the 
gradual change in the slope, one should be cautious when using 
a single slope double-/? model which might introduce a system- 
atic error in the cluster mass measurements (as also described 
in e.g. Horner 2001). The power law fit of the ICM density 
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Fig. A.6. See caption in Fig. IA.5I 



distributions at rsoo for the present sample gives an average of 

H e oc r -2.21±0.(K 

We made an additional attempt to check the difference 
in the mass modeling using a single-/? model and a double- 
ts model for CCCs caused by the steepening of the slope on 
large radial scales. We used a single-/? model with the value of 
the core radius fixed to results for the outer component in the 



double-/? model to model the electron number density distribu- 
tion. It often gives a higher value for the slope than the free 
fitting slope of the double-/? model by a few per cent, and up to 
10 per cent for the most pronounced CCC. Therefore the clus- 
ter total mass Mf m can be higher by up to 20 per cent due to the 
steepening ICM density in the outskirts. To suppress this mass 
bias, we performed the modeling as follows. We first fitted the 
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outer region (> 0.5r5oo) using a single-/? model. We then fit- 
ted the whole radial range with a single slope, double-/? model 
with its slope fixed to the best-fit value for the > 0.5r5oo region 
(cf. Appendix IA. 3b . As the slope is locally fitted, there should 
not be a significant bias due to the beta slope steepening in 
our X-ray mass modeling. With deeper observations, the thing 
that can possibly be improved in the X-ray mass modeling to 
reduce the systematics is the slope of the temperature profile 
at r5oo- However, the Monte-Carlo error estimation method ac- 
counts for this to some extent (also cf. Croston et al. 2007). 

B.6. Scaled cooling time profiles 

The cooling time is derived by the total energy of the gas di- 
vided by the energy loss rate 



where A, n g , n e and T are the radiative cooling function, gas 
number density, electron number density and temperature, re- 
spectively. We compute the upper limit of the age of the clus- 
ter as an integral from the cluster redshift z up to z = 100. 



Cooling regions are those showing cooling time less than the 
upper limit of the cluster age. The boundary radius of such a 
region is called the cooling radius. The cooling radius is zero 
when the cooling time in the cluster center is larger than the 
upper limit of the cluster age. The cooling time at the resolved 
inner most radii of the surface brightness profiles and cooling 
radii are given in Table IB. ll 

In Fig. IB. 31 we show the cooling time profiles with their 
radii scaled to r 5 oo- The CCCs show larger cooling radii in 
units of r^oo, and show much shorter central cooling times 
than the non-CCCs. The scaled cooling time profiles show a 
structurally similar behavior above 0.2rgoo for the sample. The 
best-fit power law above 0.2rsoo gives f c (r) oc r i-58±o.oi f or ^ 
whole sample. The cooling time profiles within 0.2r5oo show 
large scatter spreading over 2 magnitudes without showing 
pronounced CCC and non-CCC bimodal distribution. Similar 
studies for nearby clusters based on Chandra data can be found 
in Sanderson et al. (2006) and Dunn & Fabian (2006). 

The cooling time is calculated towards the cluster center to 
the inner most bin for which the surface brightness profile is 
resolved. There can be large uncertainties in the cooling time 
at the inner most bins where the temperature measurements are 
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Fig. A.8. Mass profiles with their mass data points binned as the temperature measurements to demonstrate the error bars. 



not resolved. The spatial resolution is important to calculate 
the proper cooling time in the cluster center. When the cluster 
core is not well resolved, the value calculated from the mea- 
sured temperature and electron number density gives the upper 
limit of the cooling time, especially for the CCCs, which show 
sharply peaked surface brightness profiles. 



B. 7. Scaled entropy profiles 

The entropy is the key to the understanding of the thermody- 
namical history since it results from shock heating of the gas 
during cluster formation. The observed entropy is generally de- 
fined as S = Tn~ 2/3 for clusters (e.g. Ponman et al. 1999). 
Radiative cooling can raise the entropy of the ICM (e.g Pearce 



Zhang et al.: LoCuSS: comparison of observed X-ray and lensing galaxy cluster scaling relations with simulations 



39 




r (Mpc) 



r (Mpc) 



Fig. A.9. See caption in Fig. lA.8 



et al. 2000) or produce a deficit below the S-T relation from show slightly lower central entropies compared to the S-T 

the self-similar prediction (e.g. Lloyd-Davies et al. 2000). scaling law. 

According to the standard self-similar model the entropy Ponman et al. (2003) suggested to scale the entropy as S oc 

scales as S oc T (e.g. Ponman et al. 1999). We investigated the T 65 based on observations for nearby clusters. We show the 

entropy-temperature relation (S-T) using So.2r 50 o as tne central S-T relations in Fig. lB.4l and Table lB.2l The Sq.2 I500 -T relation 



entropy and To.2 



0.5r 50 o 



as the cluster temperature. The CCCs (So.2r 50 o K T 



1.32±0.16 



) for our sample agrees better with the stan- 
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Fig. A.10. See caption in Fig. lA.8 



dard self-similar model S oc T. The CCCs (S 



().2i- 5l 



oc T 



0.63±0.20 



show better agreement with the empirically determined scaling 
(Ponman et al. 2003) S oc T* - 65 . The S—T relation becomes 
less steep using the entropies at OJrsoo (Fig. |B.4t . Pratt et al. 
(2006) found that further away from the cluster center, the scat- 
ter in the entropy of the S—T relation decreases. This is not ev- 
ident for our sample due to the limited radial temperature bins. 
There is no noticeable evolution in the S-T relation (e.g. us- 
ing 5 l o.3r 50 o) comparing our sample to the nearby relaxed cluster 
sample in Pratt et al. (2006) after the redshift evolution correc- 
tion due to LSS growth has been applied. 

As S oc T and M gas oc r 1 5 according to the standard self- 
similarity, one derives S oc (TM gas ) oc F°- 5 . In Fig. IB .41 we 
present both the entropy and the pressure at 0.2r 5 oo, 0.3r5oo and 
rsoo versus X-ray analog of the SZ flux. The best-fit power law 
is given in Table |B~21 The CCCs show large scatter in S—T, 
S-Yx and P-Yx relations. The P—Yx relations show least seg- 
regation between CCCs and non-CCCs, especially at rsoo- 

We scaled the radial entropy profiles using the empirically 
scaling (Ponman et al. 2003) S oc T 065 and rsoo. As shown in 
Fig. IB. 51 the scaled entropy profiles for the sample are struc- 
turally similar above 0.2r 5 oo and show least scatter around 0.2- 



0.5r5oo. After the redshift correction, we obtained an agree- 
ment between the scaled entropy profiles for this sample, the 
Birmingham-CfA sample (Ponman et al. 2003), and the cluster 
sample in Pratt et al. (2006). There is no evident evolution in 
the scaled entropy profiles after the redshift evolution correc- 
tion due to LSS growth has been applied. 

The scaled entropy profiles for our sample give the best 
fit S(r) oc r 100 ± 003 above 0.2r5oo- A similar power law was 
also found as S oc r 91 by Ettori et al. (2002), and S oc r 095 
by Piffaretti et al. (2005), respectively. The spherical accretion 
shock model predicts S oc r 1 1 (e.g. Tozzi & Norman 2001; 
Kay 2004). The combined fit of the entropy profiles for the 9 
CCCs in our sample gives S oc r 105 ± 04 i close to the prediction 
of the spherical accretion shock model and consistent with the 
relaxed nearby clusters in Pratt et al. (2006) giving S oc r 108 . 



B.8. Scaled total mass and gas mass profiles 

The mass profiles were scaled with respect to Mf m and /-500, 



respectively (Fig. IB. 61 . We found least scatter at radii above 
0.2r5oo. In the inner regions (< 0.2rgoo), the mass profiles show 
large scatter spreading one magnitude. This might be due to 
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the complex dynamics near the cluster center. As found in the 
strong lensing studies based on high quality HST data (Smith et 
al. 2005), the cluster cores can be very complicated with multi- 
clump DM halos. 

Similar to the mass profiles, the scaled gas mass profiles ap- 
pear structurally similar at radii above 0.2r5oo but show smaller 
scatter (a few per cent) than the scaled mass profiles. 

8.9. Summary of the structural similarity 

In summary, we observed the structural similarity of the scaled 
profiles of the X-ray properties outside the cluster cores. It sup- 
ports simulations that major mergers are not common recently 
(e.g. Cohn & White 2005), while the structurally similar behav- 
ior shows an empirical structural similarity, which is in agree- 
ment with simulations (e.g. Poole et al. 07) that the signatures 
of merger impact are retained by the cluster for a substantial 
fraction of cosmic time. 

There is large scatter in the cluster cores for the scaled 
profiles of the X-ray properties, but no pronounced CCC and 
non-CCC bimodal distributions. The peculiarities in the cluster 
cores for the individual clusters may explain the large scatter in 
the scaled profiles within cluster cores where the ICM does not 
always trace DM. 

As found in Majerowicz et al. (2002), the cluster total mass 
for the pronounced CCC Abell 1835 is lower by less than lOper 
cent than the mass derived ignoring the temperature drop in the 
cluster center. Therefore the cluster total masses on large radial 
scales are relatively insensitive to the cluster cores. The scaling 
relations can be less sensitive to the peculiarities in the cluster 
cores by (1) measuring the global temperatures excluding the 
cluster cores, (2) and calculating the luminosities corrected for 
the cluster cores. 
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Table B.l. Deduced properties of 37 LoCuSS galaxy clusters. Column (1): cluster name; Col. (2): electron number density of 
the inner most radial bin; Col. (3): entropy at 0.2r5oo; Col. (4): inner most radial bin of the surface brightness profile; Col. (5): 
cooling time measured at r cen ; Col. (6): cooling radius; Col. (7): r 5 oo; Cols. (8-10): gas mass, total mass and gas mass fraction 
at r 50 o, Col. (11): cluster morphology using the cluster classification in Jones & Forman (1992), denotes "single", 1 denotes 
"primary with small secondary", 2 denotes "elliptical", 3 denotes "off-center" and 4 denotes "complex". 



Name 


n e0 


50.21-500 


M:en 




'cool 


("500 


-M gas _500 


M 


/gas,500 


Mor. 




IQ- 3 cm- 3 


keV cm 2 


kpc 


Gyr 


^500 


Mpc 


10 14 M G 


10 14 M o 






RXCJ0043.4-2037 


5.8 ± 0.4 


391 ±23 


8.76 


7.8 ±0.3 


0.11 


1.08 


0.58 ± 0.06 


4.83 ± 1.43 


0.121 ±0.066 


2 


RXCJ0232.2-4420 


14.1 ±0.1 


402 ± 16 


47.14 


3.3 ± 0.2 


0.14 


1.30 


0.89 ± 0.09 


8.43 ± 2.48 


0.106 ±0.057 


1 


RXCJ0307.0-2840 


10.2 ±0.3 


359 ± 22 


12.01 


3.6 ± 0.3 


0.15 


1.14 


0.62 ± 0.06 


5.47 ± 2.01 


0.113 ±0.083 


C 


RXCJ05 16.7-5430 


2.9 ± 0.2 


474 ± 29 


17.59 


15.5 ±0.4 


0.00 


1.19 


0.78 ± 0.09 


6.40 ± 1.83 


0.121 ±0.068 


2 


RXCJ0528.9-3927 


10.6 ± 0.4 


423 ± 25 


12.87 


4.2 ± 0.3 


0.12 


1.19 


0.87 ± 0.07 


6.40 ± 1.87 


0.135 ±0.071 


3 


RXCJ0532.9-3701 


13.2 ±0.5 


555 ± 38 


8.38 


3.3 ±0.3 


0.12 


1.13 


0.57 ± 0.07 


5.40 ± 1.59 


0.106 ±0.062 





RXCJ0547.6-3152 


5.3 ±0.2 


433 ±31 


2.32 


8.3 ±0.3 


0.10 


1.21 


0.57 ± 0.07 


5.84 ± 1.73 


0.098 ± 0.058 





RXCJ0645.4-5413 


14.9 ± 1.0 


513 ±38 


1.78 


2.7 ± 0.3 


0.12 


1.26 


0.90 ±0.11 


6.60 ± 1.95 


0.136 ±0.080 


2 


RXCJ0658.5-5556 


6.1 ±0.2 


495 ± 26 


8.84 


8.2 ± 0.2 


0.11 


1.42 


1.76 ±0.17 


10.99 ±5.31 


0.160 ±0.178 


3 


RXCJ0945.4-0839 


3.8 ±0.2 


500 ± 37 


3.37 


10.3 ± 0.2 


0.05 


1.27 


0.49 ± 0.07 


6.73 ± 1.95 


0.073 ± 0.043 


2 


RXCJ0958.3-1103 


21.6 ±0.9 


337 ± 22 


2.56 


1.4 ±0.3 


0.16 


0.99 


0.43 ± 0.05 


3.29 ± 0.97 


0.132 ±0.077 


2 C 


RXCJ2129.6+0005 


50.4 ± 1.4 


332 ± 22 


2.36 


0.6 ± 0.3 


0.16 


1.06 


0.56 ± 0.06 


4.33 ± 1.27 


0.130 ±0.073 


C 


RXCJ2218.6-3853 


6.2 ± 0.3 


450 ± 32 


2.22 


6.1 ±0.3 


0.12 


1.14 


0.47 ± 0.06 


4.83 ± 1.42 


0.097 ± 0.058 


2 


RXCJ2234.5-3744 


6.0 ±0.3 


580 ± 45 


3.33 


8.2 ±0.3 


0.08 


1.35 


0.76 ±0.11 


8.11 ±2.36 


0.094 ± 0.057 


1 


RXCJ2308.3-0211 


9.2 ± 0.3 


486 ± 49 


13.27 


5.0 ± 0.3 


0.11 


1.24 


0.65 ± 0.07 


7.42 ± 2.28 


0.088 ±0.051 





RXCJ2337.6+0016 


5.7 ± 0.3 


710 ±43 


12.59 


8.2 ± 0.3 


0.05 


1.43 


0.81 ±0.07 


10.95 ± 3.22 


0.074 ± 0.039 


3 


Abell68 


6.4 ± 0.3 


487 ± 35 


3.55 


6.6 ± 0.2 


0.10 


1.21 


0.68 ± 0.07 


6.51 ± 1.93 


0.105 ±0.060 


2 


Abellll5 


35.4 ± 1.4 


430 ± 19 


2.06 


0.8 ± 0.2 


0.12 


1.07 


0.69 ± 0.04 


4.23 ± 1.14 


0.163 ±0.072 


4 C 


Abell209 


5.7 ±0.2 


369 ± 26 


3.74 


6.3 ± 0.3 


0.14 


1.15 


0.78 ± 0.08 


5.33 ± 1.71 


0.146 ±0.089 





Abell267 


7.4 ± 0.4 


364 ± 32 


3.29 


4.9 ± 0.4 


0.15 


1.06 


0.47 ± 0.04 


4.29 ± 1.30 


0.109 ±0.060 


2 


Abell383 


54.8 ± 1.5 


302 ± 14 


1.98 


0.5 ± 0.2 


0.16 


0.98 


0.33 ± 0.04 


3.17 ±0.94 


0.104 ±0.061 


C 


Abell773 


7.4 ± 0.4 


433 ± 40 


3.15 


5.7 ±0.5 


0.14 


1.33 


1.05 ±0.12 


8.30 ±2.45 


0.126 ±0.070 





Abell781 


1.8 ± 0.1 


536 ± 61 


6.28 


22.1 ± 1.1 


0.00 


1.05 


0.48 ± 0.06 


4.49 ± 1.33 


0.107 ± 0.064 


4 


AbelryOl 


1 13.6 ±1.7 


ico i 1 a 

158 ± 10 


1 -70 

1.78 


A 1 i A A 

0.2 ± 0.4 


A O A 

0.24 


A A A 

0.99 


A A A i A A-7 

0.44 ± 0.07 


i m , A A£ 

3.22 ± 0.96 


A 1 1-7 i A AOC 

0.137 ± 0.085 


A 




Abell963 


14.3 ±0.8 


393 ± 18 


2.14 


2.6 ± 0.2 


0.14 


1.14 


0.62 ± 0.06 


5.19 ± 1.52 


0.120 ±0.066 





Abelll413 


25.3 ± 1.0 


400 ± 23 


1.59 


1.7 ±0.3 


0.14 


1.18 


0.68 ± 0.08 


5.38 ± 1.57 


0.127 ±0.073 


2 


Abelll689 


30.0 ± 0.9 


505 ± 24 


1.95 


1.5 ±0.2 


0.13 


1.44 


1.05 ±0.14 


10.26 ± 2.96 


0.102 ±0.060 





Abelll758 


2.4 ± 0.2 


487 ± 39 


3.80 


16.7 ± 0.3 


0.00 


1.43 


0.85 ±0.14 


11.15 ±3.37 


0.077 ± 0.050 


2 


Abelll763 


6.5 ± 0.3 


402 ± 47 


4.00 


6.1 ±0.5 


0.13 


1.12 


0.88 ± 0.08 


4.96 ± 1.46 


0.178 ±0.094 


1 


Abelll835 


63.9 ± 1.3 


351 ±22 


2.50 


0.6 ± 0.3 


0.18 


1.30 


1.18 ±0.14 


8.01 ±2.32 


0.148 ±0.084 


C 


Abelll914 


10.5 ±0.5 


629 ± 56 


1.84 


4.7 ± 0.3 


0.11 


1.71 


1.19 + 0.13 


16.76 ±4.87 


0.071 ±0.039 


3 


Abell2204 


143.7 ± 2.3 


401 ± 28 


1.68 


0.2 ± 0.3 


0.17 


1.17 


0.84 ±0.10 


5.31 ± 1.52 


0.158 ±0.086 


C 


Abell2218 


5.8 ±0.2 


444 ± 50 


3.26 


8.0 ±0.5 


0.12 


1.07 


0.61 ±0.06 


4.18 ± 1.27 


0.147 ±0.085 





Abell2261 


12.4 ± 0.9 


408 ± 49 


3.95 


3.4 ± 0.5 


0.14 


1.19 


0.69 ±0.10 


6.02 ± 1.73 


0.114 ±0.068 





Abell2390 


40.3 ± 1.6 


484 ± 33 


2.31 


0.9 ± 0.3 


0.14 


1.29 


1.22 ±0.16 


7.67 ± 2.28 


0.159 ±0.095 


C 


Abell2667 


40.8 ± 1.4 


398 ± 27 


2.33 


0.8 ± 0.3 


0.15 


1.19 


0.77 ± 0.09 


6.02 ± 1.74 


0.128 ±0.073 


C 


Z7160 


65.1 ± 1.5 


218 ± 13 


2.53 


0.5 ± 0.3 


0.21 


0.87 


0.39 ± 0.03 


2.39 ± 0.70 


0.161 ±0.084 





Average 


















0.121 ±0.004 





1. The errors are calculated using the uncertainties of the electron number density and temperature measurements by Monte Carlo simulations. 
For each cluster, we simulated electron density and temperature profiles of a sample of 100 clusters using observed electron density and 
temperature profiles. The mass profiles and other properties of the 100 simulated clusters were calculated to derive the errors. In Zhang et al. 
(2006), the errors are simply propagated for the REFLEX-DXL clusters. 

2. c in Col. (11) denotes CCCs. 



44 Zhang et al.: LoCuSS: comparison of observed X-ray and lensing galaxy cluster scaling relations with simulations 

Table B.2. Power law, Y -Yq X 7 , parameterized entropy-temperature, entropy- Y x and pressure-F x relations. 



X 




¥ 


Y 


r 




Scatter (dex) 




Sample 












mean 
X Y 


average 
X Y 




1 0.2-0.5* 


500 


?7T^ 


jgl.55+0.13 
jgl.45+0.15 
1Q2.09+0.17 


1.32 + 0.16 
1.45 ±0.17 
0.63 ± 0.20 


0.03 
0.03 
0.31 


0.03 
0.03 
0.18 


0.05 
0.04 
0.15 


0.06 
0.06 
0.09 


whole sample 
non-CCC subsample 
CCC subsample 


keV 




keVcm 2 C W 


' 0.2-0.5j 


500 




jq2.01±0.10 
jgl.87+0.11 
jq246±0.14 


0.97 + 0.11 
1.14 ± 0.13 
0.44 ±0.17 


0.03 
0.03 
0.13 


0.03 
0.03 
0.16 


0.05 
0.04 
0.13 


0.05 
0.04 
0.06 


whole sample 
non-CCC subsample 
CCC subsample 


keV 








500 


iV 500 p^-\4/3 


jq2. 57+0.22 
1Q244+0.27 
JQ2.90+0.27 


1.00 ±0.26 
1.18 ±0.32 
0.54 ± 0.32 


0.21 
0.04 
0.38 


0.23 
0.23 
0.41 


0.10 
0.09 
0.23 


0.22 
0.22 
0.22 


whole sample 
non-CCC subsample 
CCC subsample 


keV 




keVcm 2 C, - Z "' 


v., 
M©keV 


E(z) 


4 °- 2 '500 

keVcm 2 E W 


jq— 5.53+1.27 
ig-6.69+1.61 
jq-5.63±0.90 


0.56 ± 0.09 
0.64 ±0.11 
0.22 ± 0.06 


0.07 
0.18 
0.29 


0.10 
0.03 
0.10 


0.15 
0.12 
0.45 


0.08 
0.08 
0.10 


whole sample 
non-CCC subsample 
CCC subsample 


1V10K.C V 


E(z) 


KCV tm 


jg-2.50+0.82 
jg-3.80+1.07 
jg6.25i0.79 


0.36 ± 0.06 
0.45 ± 0.07 
1.50 ±0.05 


0.09 
0.08 
0.32 


0.03 
0.03 
0.14 


0.17 
0.15 
0.42 


0.06 
0.06 
0.07 


whole sample 
non-CCC subsample 
CCC subsample 


r X 

1V10K.C V 


E(z) 


kfe £ ( z > 4/3 

keV cm 


j q — 1 . 91+1.44 
jg-3.15il.87 
jg0.67il.53 


0.36 ±0.10 
0.54 ± 0.32 
0.18±0.10 


0.09 
0.11 
0.22 


0.23 
0.23 
0.42 


0.68 
0.23 
0.27 


0.22 
0.21 
0.22 


whole sample 
non-CCC subsample 
CCC subsample 


'X 
M©keV 


E(z) 


^°- 2r 50o Fi-Y 1 
keVcm-3 C ^ 


jg-9.26il.31 
jq-8.43+1.50 
jg- 10.90+2.83 


0.50 ± 0.09 
0.44 ±0.10 
0.61 ±0.19 


0.05 
0.08 
0.22 


0.08 
0.08 
0.20 


0.18 
0.20 
0.17 


0.10 
0.10 
0.10 


whole sample 
non-CCC subsample 
CCC subsample 


Y X 


E(z) 




jQ-8.95il.02 
jQ-8.15il.18 
j g- 10.4 1 + 1. 67 


0.46 ± 0.07 
0.41 ±0.08 
0.56 ±0.11 


0.10 
0.08 
0.24 


0.08 
0.08 
0.16 


0.16 
0.17 
0.17 


0.08 
0.08 
0.09 


whole sample 
non-CCC subsample 
CCC subsample 


1V10K.C V 


kev cm 


yx 


E(z) 


r '-5oo pr-r 2 


jQ-8.32il.88 
jg-8.47+2.10 
1Q-8.47+3.45 


0.35 ±0.13 
0.36 ±0.14 
0.37 ± 0.23 


0.19 
0.25 
0.37 


0.08 
0.08 
0.27 


0.36 
0.38 
0.50 


0.13 
0.12 
0.18 


whole sample 
non-CCC subsample 
CCC subsample 


MokcV 


keVcm- 3 


r 
»500 




S FC--14/3 [ ' - 2 -»- 5 '500 \ 
keVcm 2 ~ w \^ lOfeV j 


j q3 .46+0.01 

jg3 .46+0.02 
jq3.46±0.02 


1.00 ±0.03 

0.99 ± 0.04 
1.05 ±0.04 










whole sample 

non-CCC subsample 
CCC subsample 



Here we provide both the mean scatter derived by a Gaussian fit to the histogram and the average scatter. 



Zhang et al.: LoCuSS: comparison of observed X-ray and lensing galaxy cluster scaling relations with simulations 45 




Fig. B.l. Left: Scaled radial temperature profiles for the sample. The CCCs are marked by triangles. The colors have the same 
meaning as those in Fig.Q] Right: An average temperature profile of the sample of 37 LoCuSS clusters (black, solid) compared 
to the temperature profile ranges in Markevitch et al. (1998, gray, solid), Vikhlinin et al. (2005, gray, dashed), and Pratt et al. 
(2007, gray, dash-dotted). 




0.01 0.1 1 0.2 0.5 1 



R / r 500 r / r 500 

Fig. B.2. Left: Surface brightness profile fits scaled according to the empirical scaling, Sx x T 1 '. Right: Power law slope of the 
electron number density profiles. The colors have the same meaning as those in Fig.Q] 




Fig. B.4. Entropy (upper and middle panels) and pressure (lower panels) versus temperature (upper panels) and X-ray analog of 
the SZ flux (middle and lower panels). The lines denote the best fit for our sample. The CCCs are marked by triangles. The colors 
have the same meaning as those in Fig.Q] 



Zhang et al.: LoCuSS: comparison of observed X-ray and lensing galaxy cluster scaling relations with simulations 



47 




Fig. B.5. Left: Scaled entropy profiles for the sample and the combined best fit. Right: Scaled entropy profile fits for the sample 
compared to the sample in Pratt et al. (2006, gray, hatched). The CCCs are marked by triangles in the left panel. The colors have 
the same meaning as those in Fig.Q] Tio denotes 7o.2-o.5r 5C ,o/10 keV. 




Fig. B.6. Scaled total mass profiles (left) and gas mass profiles (right). The CCCs are marked by triangles. The colors have the 
same meaning as those in Fig. [T] The error bars (a few per cent) of the gas mass profiles are too small to be seen in the right 
panel. 




Fig. B.7. Normalized cumulative cluster number count as a function of the fraction of the total bolometric luminosity attributed 
to the < 0.2r5oo region. 
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Table C.l. Deduced properties of 37 LoCuSS galaxy clusters at the radii (r 5( 5, ) determined by the method described in § [4] 
Column (1): cluster name; Col. (2): r 5( *' ; Cols. (3): X-ray mass at r 5 ^ . Cols. (4-5): weak lensing masses at r 5 ( *' for the D06 
and B07 clusters. 



Name 


Vv X 
'500 


M x (. 


500 > 




V 5 oo > 






Mpc 


10 14 M o 




10 14 
D06 




B07 


RXCJ0043. 4-2037 


1.07 


4.78 = 


= 1.40 








— 


RXCJ0232.2-4420 


1.22 


7.66 = 


= 2.20 








— 


RXCJ0307.0-2840 


1.15 


5.51 = 


= 2.01 








— 


RXCJ05 16.7-5430 


1.20 


6.49 ± 1.83 








— 


RXCJ0528.9-3927 


1.23 


6.80 = 


= 1.98 








— 


RXCJ0532.9-3701 


1.08 


5.16 = 


= 1.51 








— 


RXCJ0547.6-3152 


1.13 


5.39 = 


= 1.60 








— 


RXCJ0645.4-5413 


1.32 


6.97 = 


= 2.07 








— 


RXCJ0658.5-5556 


1.55 


14.38 = 


= 6.07 








— 


RXCJ0945.4-0839 


1.09 


5.44 = 


= 1.57 








— 


RXCJ0958.3-1103 


1.05 


3.48 = 


= 1.02 








— 


RXCJ2129.6+0005 


1.13 


4.60=; 


= 1.34 


5.51 ±3.87 




— 


RXCJ2218.6-3853 


1.09 


4.59 = 


= 1.34 








— 


RXCJ2234.5-3744 


1.26 


7.44 = 


= 2.13 








— 


RXCJ2308.3-0211 


1.19 


7.05 = 


= 2.16 








— 


RXCJ2337.6+0016 


1.18 


8.41 = 


= 2.41 








— 


Abell68 


1.16 


6.20 = 


= 1.82 


11.41 = 


= 8.49 


4.68 


± 1.25 


Abellll5 


1.31 


5.79 = 


= 1.53 


6.30 = 


= 3.99 




— 


Abell209 


1.19 


5.52 i 


= 1.74 


8.41 ±5.09 


5.08 


± 1.29 


Abell267 


0.99 


3.83 = 


= 1.16 


7.53 ±4.13 


2.27 


±0.92 


Abell383 


1.02 


3.30 = 


= 0.97 






2.95 


±0.92 


Abell773 


1.28 


7.84 = 


= 2.25 


10.01 ±5.14 




— 


Abell781 


1.10 


4.75 = 


= 1.39 


7.69 = 


= 4.27 




— 


Abell901 


1.00 


3.27 = 


= 0.96 








— 


Abell963 


1.16 


5.28 = 


= 1.52 






3.31 


±0.57 


A hell 141 3 


1 93 


5.63 i 


= 1.63 










Abelll689 


1.32 


9.24 = 


= 2.69 






11.54 


±2.13 


Abelll758 


1.18 


8.46 h 


= 2.55 


13.34 ±6.91 






Abelll763 


1.25 


6.11 = 


= 1.77 


7.64 ±3.81 


8.03 


± 1.66 


Abelll835 


1.36 


8.43 = 


= 2.45 


8.33 ±4.17 


13.41 


±2.62 


AbeM914 


1.37 


11.56 = 


= 3.30 


7.40 ± 4.23 






Abell2204 


1.29 


5.90 = 


= 1.69 


7.76 ± 5.00 






Abell2218 


1.13 


4.47 = 


= 1.36 






6.63 


± 1.36 


Abell2261 


1.22 


6.21 = 


= 1.76 


6.09 ± 3.43 






Abell2390 


1.39 


8.25 = 


= 2.41 


9.56 ± 4.39 


7.52 


± 1.56 


Abell2667 


1.24 


6.31 = 


= 1.81 










Z7160 


1.04 


2.85 = 


= 0.83 


7.29 ± 4.46 







"D06" denotes the clusters in the D06 subsample, and "B07" the B07 subsample, respectively. We used the error of AfijL as the error for the 
weak lensing masses here. 
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Table C.2. Deduced properties of 37 LoCuSS galaxy clusters at the radii (r^ w ) determined by the method described in § 16.2.31 
Column (1): cluster name; Col. (2): rj () x ' wl ; Cols. (3): X-ray mass at rj g ,wl . Cols. (4-5): weak lensing masses at rJ o' wl for the D06 
and B07 clusters. 



Name 


Yy wl 
'500 


M X( Vx.wK 
m V500 > 




M wl ( 


Vx.wk 
'500 > 






Mpc 


10 


l4 M Q 




10 14 M Q 














D06 




B07 


RXCJ0043. 4-2037 


1.06 


4.72 ± 


1.40 








— 


RXCJ0232.2-4420 


1.20 


7.54 ± 


2.20 








— 


RXCJ0307.0-2840 


1.14 


5.44 ± 


2.01 








— 


RXCJ0516.7-5430 


1.18 


6.37 ± 


1.83 








— 


RXCJ0528.9-3927 


1.22 


6.68 ± 


1.98 








— 


RXCJ0532.9-3701 


1.07 


5.10± 


1.51 








— 


RXCJ0547.6-3152 


1.12 


5.32 ± 


1.60 










RXCJ0645.4-5413 


1.30 


6.89 ± 


2.07 










RXCJ0658.5-5556 


1.53 


13.77 ± 


6.07 










RXCJ0945. 4-0839 


1.08 


5.36 ± 


1.57 










RXCJ0958.3-1103 


1.04 


3.44 ± 


1.02 










RXCJ2 129.6+0005 


1.12 


4.55 ± 


1.34 


5.45 ± 3.87 






RXCJ22 18.6-3853 


1.08 


4.53 ± 


1.34 










RXCJ2234.5-3744 


1.24 


7.34 ± 


2.13 










RXCJ2308.3-0211 


1.18 


6.96 ± 


2.16 










RXCJ2337.6+0016 


1.17 


8.28 ± 


2.41 










Abell68 


1.15 


6.13 ± 


1.82 


11.28 ± 


8.49 


4.64 


± 1.25 


Abellll5 


1.29 


5.68 ± 


1.53 


6.23 ± 


3.99 






Abell209 


1.18 


5.46 ± 


1.74 


8.33 ± 


5.09 


5.03 ± 1.29 


Abell267 


0.98 


3.77 ± 


1.16 


7.45 ± 


4.13 


2.26 ± 0.92 


Abell383 


1.01 


3.27 ± 


0.97 






2.92 


±0.92 


Abell773 


1.27 


7.72 ± 


2.25 


9.92 ± 


5.14 






Abell781 


1.09 


4.67 ± 


1.39 


7.59 ± 


4.27 






Abell901 


0.99 


3.24 ± 


0.96 










Abell963 


1.15 


5.21 ± 


1.52 






3.29 


±0.57 


Abelll413 


1.22 


5.57 ± 


1.63 










Abelll689 


1.31 


9.14 ± 


2.69 






11.43 


±2.13 


Abelll758 


1.17 


8.34 ± 


2.55 


13.19 ± 


6.91 






Abelll763 


1.24 


5.99 ± 


1.77 


7.56 ± 


3.81 


7.93 


± 1.66 


Abelll835 


1.35 


8.34 ± 


2.45 


8.25 ± 


4.17 


13.25 


±2.62 


Abelll914 


1.36 


11.35 ± 


3.30 


7.34 ± 


4.23 






Abell2204 


1.28 


5.82 ± 


1.69 


7.67 ± 


5.00 






Abell2218 


1.12 


4.41 ± 


1.36 






6.58 


± 1.36 


Abell2261 


1.21 


6.14 ± 


1.76 


6.04 ± 


3.43 






Abell2390 


1.37 


8.15 ± 


2.41 


9.46 ± 


4.39 


7.45 


± 1.56 


Abell2667 


1.23 


6.23 ± 


1.81 










Z7160 


1.03 


2.82 ± 


0.83 


7.20 ± 


4.46 







"D06" denotes the clusters in the D06 subsample, and "B07" the B07 subsample, respectively. We used the error of M™ ' as the error for the 
weak lensing masses here. 
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Table C.3. Deduced properties of 37 LoCuSS galaxy clusters at the radii (f 5 ™ sl ) determined by the method described in § 16.2.41 
Column (1): cluster name; Col. (2): r^' 81 ; Cols. (3): X-ray mass at r 5 ^ sl . Cols. (4-5): weak lensing masses at rj ( ^' sl for the D06 
and B07 clusters. 



Name 


Yv si 
'500 


, ,\- , Y v SI \ 

M V 5 oo ) 




M wl (i 


Y v sU 
500 > 






Mpc 


10 14 M o 




10 14 
D06 


M e 


B07 


RXCJ0043. 4-2037 


1.15 


5.21 ± 1.52 




— 




— 


RXCJ0232.2-4420 


1.30 


8.39 ± 2.42 




— 




— 


RXCJ0307.0-2840 


1.23 


6.02 ±2.16 




— 




— 


RXCJ0516.7-5430 


1.29 


7.24 ± 2.04 




— 




— 


RXCJ0528.9-3927 


1.33 


7.67 ±2.21 




— 




— 


RXCJ0532.9-3701 


1.15 


5.53 ± 1.62 




— 




— 


RXCJ0547.6-3152 


1.21 


5.84 ± 1.73 




— 




— 


RXCJ0645.4-5413 


1.42 


7.54 ± 2.23 




— 




— 


RXCJ0658.5-5556 


1.69 


18.47 ±6.80 




— 




— 


RXCJ0945. 4-0839 


1.17 


6.02 ± 1.73 




— 




— 


RXCJ0958.3-1103 


1.13 


3.75 ± 1.09 




— 




— 


RXCJ2 129.6+0005 


1.21 


4.97 ± 1.44 


5.89 


± 3.87 




— 


RXCJ22 18.6-3853 


1.17 


4.97 ± 1.45 




— 




— 


RXCJ2234.5-3744 


1.35 


8.11 ±2.31 




— 




— 


RXCJ2308.3-0211 


1.28 


7.65 ± 2.32 




— 




— 


RXCJ2337.6+0016 


1.26 


9.26 ± 2.67 




— 




— 


Abell68 


1.24 


6.71 ± 1.97 


12.27 


± 8.49 


4.98 


± 1.25 


Abellll5 


1.42 


6.51 ± 1.73 


6.77 


±3.99 




— 


Abell209 


1.27 


5.97 ± 1.82 


8.98 


±5.09 


5.43 


± 1.29 


Abell267 


1.06 


4.25 ± 1.27 


8.08 


±4.13 


2.40 


±0.92 


Abell383 


1.10 


3.53 ± 1.03 






3.18 


±0.92 


Abell773 


1.38 


8.76 ± 2.52 


10.71 


±5.14 




— 


Abell781 


1.19 


5.19 ± 1.53 


8.30 


±4.27 




— 


Abell901 


1.07 


3.51 ± 1.03 








— 


Abell963 


1.24 


5.75 ± 1.65 






3.47 


±0.57 


Abelll413 


1.33 


6.07 ± 1.75 










Abelll689 


1.41 


9.97 ±2.91 






12.29 


±2.13 


Abelll758 


1.26 


9.33 ±2.80 


14.35 


±6.91 






Abelll763 


1.35 


6.93 ± 1.99 


8.18 


±3.81 


8.71 


± 1.66 


Abelll835 


1.46 


9.08 ± 2.59 


8.86 


±4.17 


14.53 


±2.62 


Abelll914 


1.47 


13.04 ± 3.74 


7.87 


±4.23 






Abell2204 


1.39 


6.35 ± 1.81 


8.27 ± 5.00 






Abell2218 


1.21 


4.88 ± 1.48 






6.99 


± 1.36 


Abell2261 


1.32 


6.74 ± 1.95 


6.51 


± 3.43 






Abell2390 


1.49 


8.88 ± 2.58 


10.20 


±4.39 


7.95 


± 1.56 


Abell2667 


1.33 


6.82 ± 1.95 










Z7160 


1.12 


3.08 ± 0.88 


7.85 


±4.46 







"D06" denotes the clusters in the D06 subsample, and "B07" the B07 subsample, respectively. We used the error of M™ ' as the error for the 
weak lensing masses here. 
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C.l. See caption in Fig.|5]except that the radius (^ a' w1 ) is determined by the method described in § |6.2.3I 
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Fig. C.3. See caption in Fig.[3]except that the radius (r 50 Q' SI ) is determined by the method described in § |6.2.4| but combining the 
Yx(r) profile and the M-Yx relation from simulations in Nagai et al. (2007b). 




Fig. C.4. See caption in Fig.[6]except that the radius (/ 50 q' s1 ) is determined by the method described in § |6.2.4| but combining the 
Yx(r) profile and the M-Yx relation from simulations in Nagai et al. (2007b). 




Fig. C.5. Luminosity-lx relation. The CCCs are marked by triangles. The clusters appearing "primary with small secondary" 
or "off-center" morphology are in green, appearing "elliptical" or "complex" morphology are in blue, and appearing "single" 
morphology are in red using the classification in Jones & Forman (1992). The Fx parameters are measured at /"J q' X in the upper 

panels, at r™ wl in the left lower panel, and at rj g sl in the right lower panel. In §|H § |6.2.3| and § 16.2.41 we gave the description of 
the derivation of the radius by combining the Fx(r) profile and the M-Yx relation. 



